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ABSTRACT 

We derive a projected 2D mass map of the well studied galaxy cluster A1689 based on an entropy- 
regularized maximum-likelihood combination of the lens magnification and distortion of red back- 
ground galaxies registered in deep Subaru images. The method is not restricted to the weak regime 
but applies to the whole area outside the tangential critical curve, where non-linearity between the 
surface mass-density and the observables extends to a radius of a few arcminutes. The known strong 
lensing information is also readily incorporated in this approach, represented as a central pixel with 
a mean surface density close to the critical value. We also utilize the distortion measurements to 
locally downweight the intrinsic clustering noise, which otherwise perturbs the depletion signal. The 
resulting 2D map shows that the surface density of A1689 is smoothly varying and symmetric, similar 
to the distribution of cluster members, with no apparent substructure at r ^ 130kpc/ft. (~ 1'). The 
projected mass profile continuously steepens with radius and is well fitted by the Navarro-Frenk- White 
model, but with a surprising large concentration Cvir — I3.4I3 3, lying far from the predicted value of 
Cvir ^ 5, corresponding to the measured virial mass, M^-„ = (2.1 ± 0.2) x IQ^^Mq, posing a challenge 
to the standard assumptions defining the ACDM model. We examine the consistency of our results 
with estimates derived with the standard weak lensing estimators and by comparison with the inner 
mass profile obtained from strong lensing. All the reconstructions tested here imply a virial mass in 
the range, Mvir = (1.5 — 2.1) x lO^^M©, and the combined ACS and Subaru-2D mass reconstruction 
yields a tight constraint on the concentration parameter, Cvir = 12.7 ± 1 ± 2.8 (c2oo ~ 10), improving 
upon the statistical accuracy of our earlier ID analysis. Importantly, our best fitting profile properly 
reproduces the observed Einstein radius of 45" {zg = 1), in contrast to other weak lensing work, 
reporting lower concentration profiles, which underestimate the observed Einstein radius. 

Subject headings: cosmology: observations — dark matter — galaxies: clusters: individual (A1689) 
— gravitational lensing 



1. INTRODUCTION 

Weak gravitational lensing of background galaxies pro- 
vides a unique, direct way to study the mass distribution 
of galaxy clusters (Bartelmann & Schneider 2001) via 
the systematic shape distortion of background galaxies 
(Tyson, Wenk, & Valdes 1990; Kaiser & Squires 1993; 
Schneider & Seitz 1996; Umetsu, Tada, & Futamase 
1999) and also to a lesser extent by the magnification 
of the background (Broadhurst, Taylor, & Peacock 1995; 
Taylor et al. 1998). We have examined both these ef- 
fects in our earlier work on A1689 (Broadhurst et al. 
2005a, hereafter BG5a), where we found good consistency 
between the magnitude of the radial depletion of back- 
ground red galaxies caused by lens magnification and the 
weak lensing distortion profile of the same background 
galaxy population, which we then combined to derive an 
improved mass profile. 

A limitation of the magnification technique is the in- 
trinsic clustering of the background, which for red galax- 
ies is mainly in the form of localized groups of modest 
angular size in the background field. In principle, with 
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redshift information sharp overdensities can be isolated 
in redshift and down weighted, or with sufficiently large 
number of redshifts, the shift of the magnified luminosity 
function can be utilized independently of density fluc- 
tuations (Broadhurst et al. 1995; Zhang & Pen 2005). 
The combination of all lensing related effects is of course 
desirable, leading to the derivation of the highest preci- 
sion feasible when constructing mass maps and density 
profiles. Furthermore, because lensing effects depend on 
distance, the cosmological redshift distance relation may 
be constrained via the geometric scaling of the lensing 
signal with redshift (e.g., Taylor et al. 2007; Medezinski 
et al. 2007). 

Advances in the quality of imaging encourage a closer 
examination of the empirical effects of lensing and the de- 
velopment of more comprehensive techniques to extract 
reliable high resolution information. From space, deep 
multi-color images of massive clusters can be used to 
identify many sets of multiple images per cluster (Broad- 
hurst et al. 2005b, hereafter B05b; Gavazzi et al. 2003; 
Kneib et al. 2004; Sand et al. 2004; Smith et al. 2005; 
Bradac et al. 2006). From the ground the stable prime- 
focus wide-field cameras of Subaru and CFHT are pro- 
ducing data of sufficient quality to examine weak lensing 
distortions over a wide range of radius. More recently, 
wide-field near IR cameras, such as MOIRCS on Sub- 
aru WIRCAM on CFHT, OMEGA2000 on Calar-Alto 
and WFCAM on UKIRT, may help improve the accuracy 
of photometric redshifts for many faint galaxies. How- 
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ever, it is still the case that no set of deep high quality 
wide field images exists for any massive cluster with full 
optical-IR coverage, despite all the progress of field sur- 
veys. 

A further motivation for pursuing accurate lensing 
maps is the increased precision of model predictions 
for statistical properties of cluster-sized mass halos in 
the standard Lambda cold dark matter (ACDM) model. 
Many of the free parameters of this model now rest on 
a firm empirical foundation with relatively tight con- 
straints on the index and normalization of the power 
spectrum of density perturbation and the backgroimd 
cosmological model (e.g., Spergel et al. 2003; Tegmark 
et al. 2004; Spergel et al. 2007). In this context A^-body 
simulations have become impressively comprehensive, in 
particular the recent Millennium simulation (Springel et 
al. 2005) which simulates a huge volume of 500Mpc//i, 
and has been used to predict the mass function and evo- 
lution of nearly 100,000 group and cluster sized CDM ha- 
los. This model is tightly defined and hence amenable to 
comparisons with the real Universe, particularly for the 
case of chistcrs where baryons, which are usually omitted 
from large scale simulation, are not expected to have a 
significant impact on the shape of gravitational poten- 
tial of a cluster since the high temperature of the cluster 
gas prevents efficient cooling, and hence the majority of 
baryons simply trace the gravitational potential of the 
dominant dark matter. Accurate A^-body simulations 
based on the ACDM scenario predict a relatively shallow, 
low-concentration mass profile for massive cluster halos, 
where the logarithmic gradient flattens continuously to- 
ward the center of mass (Navarro, Frcnk, & White 1997, 
hereafter NFW) with a flatter central slope than a purely 
isothermal body interior to the inner characteristic ra- 
dius, Ts ^ 100— 200kpc//i. A useful index of the degree of 
concentration, Cvir, compares the virial radius, r^ii, to Vs 
of the NFW proflle, Cvir = ry-„/rs. This prediction for the 
CDM halo Cvir-Mvir relation has been established thor- 
oughly with high resolution simulations (e.g., Navarro 
et al. 1997; Bullock et al. 2001; Wechsler et al. 2002; 
Neto et al. 2007) with some intrinsic variation related to 
the individual assembly history of a cluster (e.g., Jing & 
Suto 2000; Tasitsiomi et al. 2004; Hennawi ct al. 2007). 
In particular, the detailed A^-body millennium simula- 
tion (Springel et al. 2005) predicts a simple relationship 
between the halo mass and concentration parameter for 
halo masses in the range of galaxy groups to massive 
clusters, as quantifled by Neto et al. (2007), who found 
that the expected median value for cluster sized halos of 
Mvir ~ 1O^^M0 to be Cvir ~ 5 at z = (c2oo ~ 4), with 
a spread of the order of A log^g Cvir = 0.1. (see Johnston 
et al. 2007 for a good summary of the state of art in halo 
concentrations based on Wechsler et al. 2006 and Neto 
et al. 2007). 

In this paper we explore further methods designed to 
achieve the maximum possible lensing precision by com- 
bining all lensing information for A1689. This cluster is 
among the most massive clusters with the largest known 
Einstein radius (^ 50"), and is one of the best studied 
clusters for lensing work (Tyson et al. 1990; Tyson & 
Fisher 1995; Taylor et al. 1998; King, Clowe, & Schnei- 
der 2002; Bardeaul et al. 2005; B05a; B05b; Oguri et 
al. 2005; Halkola, Seitz, & Pannella 2006; Bardeau et 
al. 2007; Limousin et al. 2007; Medezinski et al. 2007; 



Umetsu, Takada, & Broadhurst 2007; Okura, Umetsu, & 
Futamase 2008), located at a moderately low redshift of 
z = 0.183. In B05a we developed a "model-independent" 
method^ for reconstructing the cluster mass profile us- 
ing azimuthally-averaged weak-lensing shape distortion 
and magnification bias measurements, in the wide-field, 
Subaru images. This together with many multiple im- 
ages identified in deep Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) imaging defined a 
detailed lensing based cluster mass profile out to the clus- 
ter virial radius [r ^ 2h~^ Mpc). The combined strong 
and weak lensing mass profile is well fitted by an NFW 
profile (Navarro et al. 1997) with high concentration 
of Cvir ~ 13.7, which is significantly larger than the- 
oretically expected (cVir 5) for the standard ACDM 
model (Bullock et al. 2001; Neto et al. 2007), although 
the degree of concentration is still controversial (B05a; 
Medezinski et al. 2007; Limousin ct al. 2007). Such a 
high concentration is also seen in other massive clusters 
from careful lensing work, such as MS 2137-23 (c2oo — 12, 
Gavazzi et al. 2003) and CL0024-I-1654 (caoo ^ 22, Kneib 
et al. 2003). These results could raise serious ques- 
tions regarding the basic assumptions behind the ACDM 
model. If clusters collapse earlier than predicted then it 
is expected that denser and hence more concentrated ha- 
los will develop in the context of CDM (Wechsler et al. 
2002). On the other hand, it has been argued that part 
of this discrepancy from lensing observations could be 
reconciled by observational effects such as triaxiality of 
CDM halos (Oguri et al. 2005; Hennawi et al. 2007; 
Sereno 2007; Corless & King 2007), and the projection 
of structure along the line of sight (e.g.. King & Corless 
2007), both of which boost the projected surface mass 
density and hence the lensing signal. Such observational 
biases in the lensing-based concentration parameter have 
been explored in details by Hennawi et al. (2007) on 
the basis of A^-body simulations, indicating a positive 
bias of ~ 30% in the halo concentration derived from 2D 
lensing measurements. Although A1689 is a very round 
shaped cluster with evidence of only modest substruc- 
ture (Teague, Carter, & Grey 1990; Girardi et al. 1997; 
Andersson & Madejski 2004; Czoske 2004; B05b), such a 
chance alignment of structure could be a potential source 
of high concentrations. Furthermore, for a reliable mea- 
surement of the cluster mass profile, systematic errors 
inherent in the lensing measurements, such as the un- 
certainty in the background redshift distribution and the 
dilution effect on the lensing signal due to contamination 
by cluster members (B05a; Medezinski et al. 2007), need 
to be taken into account. 

The paper is organized as follows. We briefly summa- 
rize in §2 the basis of cluster weak lensing. In §3 we de- 
scribe the observational data and the background sample 
selection for the weak lensing analysis; we then summa- 
rize our joint weak lensing analysis of shape distortion 
and magnification bias data. In §4 we present a method 
for reconstructing the two-dimensional mass distribution 
of A1689 from combined weak lensing shape distortion 
and magnification bias measurements. In §5 we derive 

^ Wc remind the reader that model dependence is unavoidable 
to some extent in scientific analysis. In this work we define the 
term "model independent" to refer to those methods without prior 
assumptions about the functional form of the lensing profiles and 
distributions. 
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mass profiles of A1689 from weak Icnsing data using 
three different methods, and compare resulting mass pro- 
files; we also combine our weak lensing mass profiles with 
strong lensing constraints from previous studies to test 
the CDM paradigm; then, we assess carefully various 
sources of potential systematic error in the halo concen- 
tration parameter derived from the lensing observations. 
Finally, summary and discussions are given in §6. 

Throughout this paper, we use the AB magnitude sys- 
tem, and adopt a concordance ACDM cosmology with 
(ilmo = 0.3, JIao = 0.7, h = 0.7). In this cosmology one 
arcminute corresponds to the physical scale 129kpc//i 
for this cluster. The reference center of our analysis is 
fixed at the center of the cD galaxy: RA = 13 : 11 : 
29.52, Dec = -01 : 20 : 27.59 (J2000.0). 

2. CLUSTER WEAK LENSING 

Weak gravitational lensing is responsible for the weak 
shape-distortion and magnification of the images of back- 
ground sources due to the gravitational field of inter- 
vening foreground clusters of galaxies and large scale 
structures in the universe. The deformation of the im- 
age can be described by the 2x2 Jacobian matrix Aaff 
{a,P = 1, 2) of the lens mapping. The Jacobian Aa0 is 
real and symmetric, so that it can be decomposed as 



■AaP = (1 - K)Sal3 - Tc/S, 



72 -7i 



(1) 
(2) 



where 5ai3 is Kronecker's delta, is the trace-free, sym- 
metric shear matrix with 7a being the components of 
spin-2 complex gravitational shear 7 7i + 472, de- 
scribing the anisotropic shape distortion, and k is the 
lensing convergence responsible for the trace-part of the 
Jacobian matrix, describing the isotropic area distortion. 
In the weak lensing limit where K,\'-f\ <C 1, Fq,^ induces 
a quadrupole anisotropy of the background image, which 
can be observed from ellipticities of background galaxy 
images. The flux magnification due to gravitational lens- 
ing is given by the inverse Jacobian determinant. 



1 



det^ (1 - k)2 - |7|2 



(3) 



where we assume subcritical lensing, i.e., detA{d) > 0. 
The lensing convergence is expressed as a line-of-sight 

projection of the matter density contrast out to the 
source plane (S) weighted by certain combination g of 
CO- moving angular diameter distances (e.g., Jain et al. 
2000), 
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where o is the cosmic scale factor, x is the co-moving 
distance; is the surface mass density of matter, 

Tim = Jo^d,x a{pm — p), with respect to the cosmic mean 
density p, and S^rit is the critical surface mass density 
for gravitational lensing. 
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with Dg, Dd, and Dds being the angular diameter dis- 
tances from the observer to the source, from the observer 



to the deflecting lens, and from the lens to the source, re- 
spectively. For a fixed background cosmology and a lens 
redshift Zd, Scrit is a function of background source red- 
shift Zg. For a given mass distribution S(0), the lensing 
signal is proportional to the angular diameter distance 
ratio, Dds/Dg. 

In the present weak lensing study we aim to reconstruct 
the dimensionless surface mass density k from weak lens- 
ing distortion and magnification data. To do this, we uti- 
lize the relation between the gradients of k and 7 (Kaiser 
1995; Crittenden et al. 2002), 



Ak(6/) = d'^d'^Tc^pie) = 2V*-f{e) 



(6) 



where T) is the complex differential operator 2? = 
{&l — 9|)/2 -I- idid2- The Green's fimction for the two- 
dimensional Poisson equation is A~^(0, 0') = ln|0 — 
0'\/{2'k), so that equation (6) can be solved to yield the 
following non-local relation between k and 7 (Kaiser & 
Squires 1993): 

K{e) = ^J (fo' D*{e - 9')j{e') (7) 

where D{d) is the complex kernel defined as 



D{0) 



(8) 



Similarly, the spin-2 shear field can be expressed in terms 
of the lensing convergence as 



7(0) 



(9) 



Note that adding a constant mass sheet to k in equa- 
tion (9) does not change the shear field 7(6) which 
is observable in the weak lensing limit, leading to the 
so-called mass-sheet degeneracy based solely on shape- 
distortion measurements (e.g., Bartelmann & Schneider 
2001; Umetsu et al. 1999). In general, the observable 
quantity is not the gravitational shear 7 but the reduced 
shear, 

in the subcritical regime where det^ > (or 1/g* in the 
negative parity region with det^ < 0). We see that the 
reduced shear g is invariant under the following global 
transformation: 

K{e) Xk{0) + 1 - a, 7(6/) ^ A7(0) (11) 

with an arbitrary scalar constant A 7^ (Schneider & 
Seitz 1995). This transformation is equivalent to scaling 
the Jacobian matrix A{d) with A, A{9) XA{d). This 
mass-sheet degeneracy can be unambiguously broken by 
measuring the magnification effects, because the magnifi- 
cation p transforms under the invariance transformation 
(11) as 

p{e) ^ X^p{d). (12) 

3. DATA ANALYSIS 

In this section we present a full technical description 
of our weak lensing distortion and magnification analyses 
on A1689 based on the Subaru images, which were anal- 
ysed in our earlier work of B05a and Medezinski et al. 
(2007). This work is also based on the same weak lensing 
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shape and magnification measurements as used in B05a. 
We note that Mcdezinski ct al. (2007) used a slightly 
different analysis pipeline so as to include the shape mea- 
surements of bright cluster members and noisier objects 
(as well as better resolved red background galaxies), op- 
timizing for the weak lensing dilution analysis including 
the measurements of cluster light and cluster luminosity 
function where the completeness is crucial; Mcdezinski et 
al. (2007) included a blue background population in the 
weak lensing shape analysis, in addition to the red back- 
ground population which B05a and this work are based 
on; The magnification information, on the other hand, 
was not taken into account in Mcdezinski ct al. (2007). 

In B05a we simply assumed that the mean redshift of 
the red background galaxies is Zg = 1. In the present 
work we improve the accuracy of determination of the 
cluster mass and concentration parameters, by taking 
into account the redshift distribution of red background 
galaxies examined by Mcdezinski et al. (2007) based on 
the multicolor photometry of Capak et al. (2004) in the 
HDF-N (see §3.3). 

3.1. Subaru Data, and Photom,etry 

For our weak lensing analysis we used 
Subaru/Suprimc-Caui imaging data of A1689 in V 
and SDSS i' retrieved from the Subaru archive, SMOKA 
(see B05a and Medezinski et al. 2007 for more details). 
The FWHM in the co-added mosaic image is 0"82 in 
V and 0"88 in i' with 0'f202 pix^i, covering a field of 
« 30' X 25'. 

Photometry is based on a combined V image with 
SExtractor (Bertin & Arnouts 1996), where the i' im- 
age is used as the source detection image. We adopt 
the following key configuration parameters in SExtrac- 
tor: DETECT _MIN AREA = 5, DETECT .THRESH = 
ANALYSIS. THRESH = 3. The hmiting magnitudes 
are obtained as V = 26.5 and i' = 26.0 for a 3ct de- 
tection within a 2"aperture. A careful background se- 
lection is critical for a weak lensing analysis so that un- 
lensed chistcr members and foreground galaxies to not di- 
lute the true lensing signal of the background (see B05a 
and Medezinski et al. 2007). We identify an E/SO se- 
quence of cluster galaxies in the color-magnitude (CM) 
diagram, which can be defined by the linear CM rela- 
tion: {V - i')E/s = -0.02094i'+ 1.255 (B05a; see Figure 
1 of Medezinski et al. 2007 for the CM diagram). For 
the number counts to measure lensing magnification, we 
define a sample of galaxies that are redder than the clus- 
ter sequence by {V — i') > 1.0 and 20 < z' < 25.5, which 
yields a total of N^j^ = 8, 907 galaxies, or the mean surface 
number density of n^u = 12.0 arcmin"^. For the magni- 
fication bias analysis, a conservative magnitude limit of 
i' < 25.5 is adopted to avoid incompleteness. 

3.2. Weak Lensing Distortion Analysis 

We use the IMCAT package developed by N. Kaiser ^ 
to perform object detection, photometry and shape mea- 
surements, following the formalism outlined in Kaiser, 
Squires, & Broadhurst (1995; hereafter KSB). We have 
modified the method somewhat following the procedures 
described in Erben et al. (2001). We used the same 

® http://www.ifa.hawaii/kaiser/IMCAT 



analysis pipeline as in B05a, Umetsu et al. (2007), and 
Okabe & Umetsu (2008). 

Object Detection 

Objects are first detected as local peaks in the image 
by using the IMCAT hierarchical peak-finding algorithm 
hfindpeaks which for each object yields object parame- 
ters such as a peak position, an estimate of the object 
size (rg), the significance of the peak detection (v). The 
local sky level and its gradient arc measured around each 
object from the mode of pixel values on a circular annulus 
defined by inner and outer radii of 16 x Vg and 32 x rg (see 
Clowe et al. 2000). In order to avoid contamination in 
the background estimation by bright neighboring stars 
and/or foreground galaxies, all pixels within 3 x rg of 
another object are excluded from the mode calculation. 
Total fluxes and half-light radii (rh) are then measm-ed 
on sky-subtracted images using a circular aperture of ra- 
dius 3-\/2 X fg from the object center. Any pixels within 
2.5 X Tg of another object are excluded from the aper- 
ture. The aperture i'-magnitude is then calculated from 
the measured total flux and a zero-point magnitude. Any 
objects with positional differences between the peak lo- 
cation and the weighted-centroid greater than d = 0.4 
pixels are excluded from the catalog. 

Finally, bad objects such as spikes, saturated stars, and 
noisy detections must be removed from the weak lensing 
object catalog. We removed from our object catalog (1) 
extremely large objects with rg > 10 pixels, (2) objects 
with low detection signiflcance, < 7, (3) objects with 
large raw ellipticities, |e| > 0.5, (4) noisy detections with 
unphysical negative fluxes, and (5) objects containing 
more than 10 bad pixels, nbad > 10. This selection pro- 
cedure yields an object catalog with N = 62, 384 (82.6 
arcmin~2). 

Weak Lensing Distortion Measurements 

To obtain an estimate of the reduced shear, ga = 
7a/(l ^ '*) {a = 1,2), we measure the image ellipticity 
Ga = {Qii - Q22, Q12} /{Qii + Q22) from the weighted 
quadrupole moments of the surface brightness of individ- 
ual galaxies defined in the above catalog, 

ga/3 = Jd^ewie)e„epi{9) (a,/? ==1,2) (13) 

where I{d) is the surface brightness distribution of an 
object, W{6) is a Gaussian window function matched to 
the size of the object. 

Firstly the PSF anisotropy needs to be corrected using 
the star images as references: 

e'a = ea- P^iq} (14) 
where Psm is the smear polarizaMlity tensor being close to 
diagonal, and g* = {Psm)a^^* stellar anisotropy 

kernel. We select bright, unsaturated foreground stars 
of 20 ^ i' ^ 22.5 identified in a branch of the half-light 
radius (r/J vs. magnitude (i') diagram to measure (7*. 
In order to obtain a smooth map of (7* which is used 
in equation (14), we divided the co-added mosaic im- 
age of 9K X 7.4K pixels into 5x4 blocks, each with 
1.8K X 1.85K pixels. The block length is based on the 
typical coherent scale of PSF anisotropy patterns. In 
this way the PSF anisotropy in individual blocks can be 
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well described by fairly low-order polynomials. We then 
fitted the q* in each block independently with second- 
order bi-polynomials, (9), in conjunction with iterative 
cr-clipping rejection on each component of the residual: 
Se* = e* - (P;„J"'3q*(0). The final stellar sample con- 
sists of 540 stars (i.e., N^, ~ 30 stars per block), or the 
mean surface number density of n* = 0.72 arcmin"^. It 
is worth noting that the mean stellar ellipticity before 
correction is (ei*,e2*) ~ (—0.013,-0.018) over the data 
field, while the residual e* after correction is reduced to 
ef«^ = (0.47 ±1.32) x 10"^, e^''^" = (0.54 ±0.94) x lO"'^; 
The mean offset from the mill expectation is reduced 
down to |e"''"| = (0.71 ± 1.12) x lO""*. On th e other 
hand, the rms value of stellar ellipticities, a,,^ = (|e*p), 
is reduced from 2.64% to 0.38% when applying the 
anisotropic PSF correction. We show in Figure 1 the 
quadrupole PSF anisotropy field as measured from stel- 
lar ellipticities before and after the anisotropic PSF cor- 
rection. Figure 2 shows the distribution of stellar ellip- 
ticity components before and after the PSF anisotropy 
correction. From the rest of the object catalog, we se- 
lect objects with fh* l^^h pixels as an i'-selected 
weak lensing galaxy sample, where fh* ~ 2.4 pixels is 
the median value of stellar half-light radii, correspond- 
ing to half the mcxiian width of circularized PSF over the 



data field. An apparent magnitude cutoff of 20 



26 



is also made to remove from the weak lensing galaxy 
sample bright foreground/cluster galaxies and very faint 
galaxies with noisy shape measurements. 

Second, we need to correct image ellipticities for the 
isotropic smearing effect caused by atmospheric seeing 
and the window function used for the shape measure- 
ments. The pre-seeing reduced shear Qa can be estimated 
from 

go. = {Pg-'Ue'^ (15) 

with the pre-seeing shear polarizability tensor -Pf^y. We 
follow the procedure described in Erben et al. (2000) 
to measure P^. We adopt the scalar correction scheme, 
namely 

{PgU = l^AP<,]Sa0 = (16) 

(Erben et al. 2000: Hoekstra et al. 1998; Hudson et al. 
1998; Okabe & Umetsu 2008). 

In order to suppress artificial effects due to the noisy 
Pg estimated for individual galaxies, we apply filtering 
to raw Pg measurements. Firstly, we discard those noisy 
objects which have negative raw Pg values. Secondly, we 
compute for each object a median value of P| among A'^- 
neighbors in the size Vg and magnitude i' plane to define 
object parameter space: for each object, A^-neighbors are 
identified in the size (r^) and magnitude {i') plane; the 
median value of Pg is used as the smoothed P^ for the 
object, {Pg), and the variance cr"^ of g = gi +ig2 is calcu- 
lated using equation (15). The dispersion Ug is used as 
an rms error of the shear estimate for individual galax- 
ies. We take A^ = 30. After filtering noisy Pg measure- 
ments, the minimum value of {Pg) is « 0.035. Figure 3 
shows the averaged {Pg) as a function of object size Vg for 
the I'-selected weak lensing galaxy sample. The mean of 
{Pg) over all galaxies in the sample is obtained as 0.307, 
mostly weighted by galaxies with = 2 — 3 pixels. The 



mean variance cr^ over the galaxy sample is obtained as 
~ 0.152, or ^Ja^ « 0.39. Finally, we use the following 

estimator for the reduced shear: = e^/ {P^)- The 
final i'-selected galaxy sample contains 30, 369 galaxies 
or fig ~ 40.3 arcmin^^. 

3.3. Red Background Selection 

As demonstrated by B05a and Medezinski et al. 

(2007), it is crucial in the weak lensing analysis to make 
a secure selection of background galaxies in order to min- 
imize contamination by cluster/foreground galaxies and 
hence to make an accurate determination of the cluster 
mass, otherwise dilution of the distortion signal results 
from the inclusion of unlensed cluster galaxies, particu- 
larly at small radius where the cluster is relatively dense. 
This dilution effect is simply to reduce the strength of 
the lensing signal when averaged over a local ensemble 
of galaxies, in proportion to the fraction of unlensed clus- 
ter/foreground galaxies whose orientations are randomly 
distributed, thus diluting the lensing signal relative to the 
reference background level derived from the background 
population (Medezinski et al. 2007). With a pure red 
background sample (B) as a reference, one can quantify 
the degree of dilution for a galaxy sample (G) contain- 
ing A"cL cluster galaxies and A'bg background galaxies 
in terms of the strengths of the averaged tangential shear 
signal {g+{0)) as (Medezinski et al. 2007) 

{gf\e)) {Das/D,)[ 



BG i- 
A^BG 



CL 



^(G) 



{g'^'m {Das/D,,,^^,^ 

(17) 

where {DdslP>s)z,>zd^ are averaged distance ratios for 
respective background populations; if the two samples 

contain the same background population, then 5d = 

{g\. )/{g+ ) — 1- The degree of dilution thus varies de- 
pending on the radius from the cluster center, increas- 
ing towards the cluster center. Medezinski et al. (2007) 
found for their green galaxy sample {[V — ^'Ie/sq-os) 
containing the cluster sequence galaxies in A1689 that 
the fraction of cluster membership, A^cl/(A^bg + -^Cl), 
tends ~ 100% within d <2'. 

For our weak k;nsing distortion analysis we define a 
sample of red background galaxies whose colors are redder 
due to large fc-corrections than the CM relation, or red 
sequence, of cluster member galaxies. These red back- 
ground galaxies are largely composed of early to mid-type 
galaxies at moderate redshifts (Medezinski et al. 2007). 
Cluster member galaxies are not expected to extend to 
these colors in any significant numbers because the in- 
trinsically reddest class of cluster galaxies, i.e. E/SO 
galaxies, are defined by the red sequence and lie blueward 
of chosen sample limit, so that even large photometric er- 
rors will not carry them into our red sample. This can 
be demonstrated readily, as shown in Figure 4, where we 
plot the mean tangential shear strength {g+), averaged 
over a wide radial range of 1' < 6* < 18', as a function of 
color limit by changing the lower color limit progressively 
blueward. Here we do not apply area weighting to en- 
hance the effect of dilution in the central region. We take 
the lower (bluer) color limit of +0.22 mag where the clus- 
ter contribution is negligible, defining the red background 
sample by A(y - i') = {V - i') - {V - i')E/so > 0.22, 
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as adopted by B05a. Figure 4 shows a sharp drop in the 
lensing signal at A{V — i') ^ 0.1, when the cluster red 
sequence starts to contribute significantly, thereby re- 
ducing the mean lensing signal. At A{V — i') ^ 0.1, the 
mean lensing signal of the red background stays fairly 
constant, ((?+) ~ 0.143, ensuring that our weak lensing 
measurements are not sensitive to this particular choice 
of the color limit (sec §5.5). Similarly, on the left of Fig- 
ure 4, we show the mean distortion strength for our blue 
galaxy sample with the magnitude limit in the interval 
23 < i' < 25.5, so as to take only faint blue galaxies. The 
lower color limit is set to /S.{V — i') > —1.5. For galaxies 
with colors bluer than the cluster sequence, cluster galax- 
ies are present along with background galaxies, since the 
cluster population extends to bluer colors of the later 
type members. Consequently, the mean lensing signal of 
the blue sample is systematically lower than that of the 
red sample, unless we take the upper (redder) color limit 
around A{V — i') ~ 0.7 where, however, the distortion 
measurement is quite noisy (Figure 4). We therefore ex- 
clude blue galaxies from our weak lensing analysis. Note 
that the background populations do not need to be com- 
plete in any sense but should simply be well defined and 
contain only background. This color-magnitude selection 
criteria yielded a total of Ng = 5728 galaxies, or the mean 
surface number density of fig = 7.6 galaxies arcmin"^. 
For the red background sample, we found C7g ~ 0.133, 

or Ri 0.36, which is slightly smaller than the rms 

dispersion for the i'-sclcctcd galaxy sample. 

We need to estimate the depths of our color-magnitude 
selected red samples when measuring the cluster mass 
profile, because the lensing signal depends on the source 
redshifts in proportion to Dds/Dg. Medezinski et al. 
(2007) utilized the multicolor photometry of Capak et 
al. (2004) based on Subaru U BVRIZ imaging covering 
0.2 dcg^, and estimated photometric redshifts for color- 
magnitude selected background galaxy samples. Using 
the Capak et al.'s photometric redshift distributions, 
Medezinski et al. (2007) estimated a mean redshift for 
the red background with i' > 18 mag to be Zs « 0.87 
at fainter magnitude limits of i' = 26-27, and also cal- 
culated weighted mean lensing depths {D^s/Ds) for re- 
spective background populations as a function of appar- 
ent i' limiting magnitude, and found that the averaged 
distance ratio {Dds/Ds) grows only slowly with increas- 
ing apparent magnitude limit: {Dds/Ds) = 0.693 ± 0.02 
for the red sample with i'^^^ = 25.5-26.5 (see Figure 9 of 
Medezinski et al. 2007), and the corresponding redshift 
equivalent to this mean distance is Zs^d = 0.68 ± 0.05, 
where the distance-equivalent redshift Zg^o is defined by 
the following equation: 

Dds \ _ Dds 

'd7/^ ~ 'd7 



(18) 



Therefore, we can safely assume that our two different 

red samples for the distortion and magnification mea- 
surements (with z^jji^ = 26.0 and 25.5, respectively) have 
nearly the same depths. We note that B05a assiimcd 
that all of the background galaxies are located at a sin- 
gle redshift of 2:^ = 1, corresponding to Dds/Dg = 0.772, 
which underestimates Scrit by ^ 11% and hence underes- 
timates the cluster mass accordingly. We will come back 
to this issue in §5.5. 



4. WEAK LENSING MAP-MAKING 

This section describes a maximum entropy method 
(MEM) for reconstructing the two-dimensional cluster 
mass distribution from combined shape distortion and 
magnification bias observations. 

As described above (see §3.3), only the red selected 
background is used for the measurement of the lens dis- 
tortion and magnification, in order to minimize the effect 
of dilution. 

4.1. Shape Distortion Data 

For map-making, we pixelize the distortion data into 

a regular grid of iVdata ~ 21 x 17 = 357 independent 
pixels covering a field of « 30' x 24'. The pixel size is set 
to A 



pix 



I'A, and the mean galaxy counts per pixel is 
figAfipix ~ 15 with AQpix = Apj^ being the solid angle 
per pixel. 

We compute an estimate ga{(^i) for the reduced shear 
9a{di) on a regular grid of cells (z 



l,2,...,iVpix) as 



/ ^ Ukga.kl ^ Uk = ga,i, (19) 
/cGcelli fcfEcelH 

where ga,k is a noisy estimate of the ath component of 
the reduced shear for the kth galaxy, and Uk = ^/<^1^k i^ 
its inverse- variance weight (see §3.2). The error covari- 
ance matrix for the pixelized reduced shear g = g-i+ ig2 
(19) is then diagonal, and given by 



1 



{Aga,iAgf3j) = -(Jg^iSaokj, 



(20) 



where tr? • is the error variance for the jth pixel defined 



as 



(21) 



Here we have used {Ag„^hAg/3^i) = (1/2)ct^ j.(5a/34; for 
individual shear estimates ga,k- The per-pixel rms dis- 
persion for g{6) is then reduced down to '^Ug/N 

0.36/VT5 - 0.092. 

In Figure 5 we show the reduced-shear field obtained 
from the red galaxy sample, where for visualization pur- 
poses the ga{d) is resampled on to a finer grid and 
smoothed with a Gaussian with FWHM = 2'. A co- 
herent tangential shear pattern is clearly seen in Figure 
5 around the cluster center. 

4.2. Magnification Bias Data 

Lensing magnification, ij,{6), influences the observed 
surface density of background sources, expanding the 
area of sky, and enhancing the observed flux of back- 
ground sources. In the subcritical regime (Broadhurst, 
Taylor & Peacock 1995; Umetsu et al. 1999), the mag- 
nification /i(0) is given by equation (3). The count-in- 
cell statistic is measured from the flux-limited red galaxy 
sample (see §3.1) on the same grid as the distortion data: 



N{< mcut; 0i) 



/cGcelli 



1 = A^, 



(22) 



with TOcut being the magnitude cutoff corresponding to 
the fiux-limit (i^ut ~ 25.5). The normalization and slope 
of the unlensed number counts A^o(< Ticut) for our red 
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galaxy sample are reliably estimated as n^^^o = 12.6 ± 
0.23 arcmin-2 and s = d\og No{m)/dm ='o.22±0.03 
from the outer region > 10' (B05a). The mean galaxy 



counts per pixel is thus A^o 



pix 



25. The 



magnification bias at the ith cell is then estimated as 1 + 
= %/ny^,o = Ni/No, with n^,i = A^^,j/Aripix, where 
the dilution effect 6d is negligible for our red-background 
sample; otherwise Ni/N^ = \ + 5^-\-5d- The slope is less 
than the lensing invariant slope, s = 0.4, and hence a net 
deficit of background galaxies is expected: 



(23) 



(Broadhurst et al. 1995; BG5a). In the limit of weak lens- 
ing where k, |7| < 1, (5^(0) ~ (5s - 2)k(0) ~ -0.9k(6') 
with s = 0.22. 

The masking effect due to bright cluster galaxies and 
bright foreground objects is properly taken into account 
and corrected for (B05a). Wc conservatively account for 
the masking of observed sky by excluding a generous area 
Tiab around each masking object, where a and h are de- 
fined as I'mask = 3 times the major (A_IMAGE) and minor 
axes (B_IMAGE) computed from SExtractor, correspond- 
ing roughly to the isophotal detection limit in our config- 
uration (see §3.1). Wc calculate the correction factor for 
this masking effect as a function of radius from the cluster 
center, and renormalize the number density of each cell 
accordingly. The masking area is negligible at large radii, 
and increases up to ^ 20% of the sky close to the clus- 
ter center, 9 ^ 3'. B05a showed that the magnification 
bias measurements with and without the masking cor- 
rection are roughly consistent with each other, and with 
the NFW prediction from the ACS strong lensing ob- 
servations (B05b) and the Subaru distortion profile (see 
Figure 2 of B05a), even though all the measurements 
have different systematics. Note that if we use the mask- 
ing factor fmask of 2 or 4, instead of 3, the results shown 
below remain almost unchanged (see §5.5 for more de- 
tails) . 

Figure 6 shows the resulting magnification-bias distri- 
bution derived from the red galaxy sample based on the 
SExtractor photometry (§3.1). A clear depletion of the 
red galaxy counts is visible in the central, high-density 
region of the cluster. On the other hand, it has been 
argued that estimates of the lensing magnification based 
on number counts suffer from noise arising from the in- 
trinsic clustering of the source galaxies (e.g., Zhang & 
Pen 2006). Indeed, some variance is apparent in the spa- 
tial distribution of red galaxies. In particular, a local 
enhancement of red background galaxies can be seen in 
Figure 6 around ^ (0',5'), having an overdensity of 
nil/ no ~ 1.5. This may explain the discrepant point at 
^ ~ 5' in the magnification bias profile of B05a. We will 
come back to this issue later in §4.5. 

4.3. Maximum Entropy Mass Reconstruction Method 

The relation between distortion and convergence is 
non-local, and the convergence k derived from distortion 
data alone suffers from a mass sheet degeneracy (see §2). 
However, by combining the distortion and magnification 
measurements the convergence can be obtained unam- 
biguously with the correct mass normalization. 

Here we combine pixelized distortion and depletion 
data of red background galaxies in a maximum likeli- 



hood sense, and reconstruct the two-dimensional distri- 
bution of the lensing convergence field, k{6). Several au- 
thors have proposed maximum likelihood and maximum 
entropy methods for reconstructing the projected mass 
distribution from joint weak lensing observations of the 
shape distortion and magnification bias effects (Bartel- 
mann 1995; Seitz, Schneider, & Bartelmann 1998; Bridle 
et al. 1998; Umetsu et al. 2007). In the present study, we 
utilize a maximum entropy method extended to account 
for positive/negative distributions of the underlying field 
starting with the method proposed in the context of 
interferometric observations of the Cosmic Microwave 
Background (CMB) radiation by Maisinger, Hobson, & 
Lasenby (1997) and Hobson & Lasenby (1998). 

Bayes' theorem states that given a hypothesis H and 
some data D the posterior probability Y'r{H\D) is the 
product of of the likelihood Vt{D\H) and the prior prob- 
ability Pr(ff), with a proper normalization by the evi- 
dence Pr(£>) (Maisinger et al. 1997): 

In the context of map-making, the hypothesis H is taken 
as the pixehzed image p = {pi} (i = 1, 2, A^pix)- Since 
the evidence in Bayes' theorem is a normalization con- 
stant for a given dataset D, we maximize Pr(p|_D) oc 
Pr(p) Pr(Z)|p) with respect to p. We assume that the er- 
rors on the data follow a Gaussian distribution, so that 
the likelihood is given as Pr(Z)|p) oc exp [— x^(p)/2] , 
with being the standard misfit statistic. Then, for a 
given form of the entropy function 5(p, m), the entropic 
prior is written as 

Pr(p) a cxp [a'5(p, m)] , (25) 

where m = {mj} {i = 1, 2, A^pix) is a set of model 
parameters for the pixelized image, p. We follow the 
prescription given in Maisinger et al. (1997) and Hobson 
& Lasenby (1998), and define the cross-entropy function 
for p as 



5(p,m) = ^ 



i=l 



ipi - 2mi - Pi In 



jij +Pi 
2mi 



(26) 



where ipi = -\/pf-i-4mf. We take pi = Ki = K{6i) as 
the image to be reconstructed, and express a set of dis- 
cretized K-values as p = {ki} [i = 1, 2, A^pix)- In gen- 
eral, an entropy rcgularization helps to reduce the sensi- 
tivity of the least solutions to small-scale noise in data, 
by imposing smoothness constraints on the solutions. 
This MEM prior ensures that k(0) — > in the noise- 
dominated regime, or low-density regions in the outskirts 
of the cluster. Note that unlike conventional MEM func- 
tions, this MEM prior is free from the "positive bias" in 
the reconstructed image (= signal -|- residual noise), and 
this bias is more significant in the low signal-to-noise 
ratio, or low density regions (r — > rvir) where we are 
interested in measuring k. 

We take into account the non-linear, but subcritical, 
regime of the lensing properties, k and 7^, for a MEM 
mass reconstruction (see Bridle et al. 1998). We use 
equations (7) and (9) to relate the gravitational shear 
and the convergence fields, 7 and k. The log posterior 
probability function, f (p) = — lnPr(p|£'), is then ex- 
pressed as a linear sum of the shear /magnification data 
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log-likelihoods (Schneider, King, & Erben 2000) and the 
entropy term (Maisinger et al. 1997): 



F{p) = -x\p)-aS{p,in), 

X'(p)=xg(p)+X'(p) 

ATpix 2 



i—l a—1 



9a,i - 9aAP) 
-I, 



X/j, — ^ ] 



[Ni-N,ip)f 



(27) 
(28) 

(29) 
(30) 



where a(> 0) is the dimensionless regularization con- 
stant, and (Tg^^ ^ is the pcT-componont rms error for the 
pixclizcd distortion measurement (a ~ 1,2), ct^j ^ = 
^92, i = cr~g,i/V^'-, ga,i{p) (equation [10]) and Ni{p) = 
Mi(p)^'^*~^-^o are the theoretical expectations for g^^i 
and Ni, respectively, where Afipix,i is the effective ob- 
served area of the ith pixel excluding the masking area 
by the cluster members. In defining x^, we have used 
the Gaussian approximation for the Poisson distribution 
of count-in-cell statistics (A^o ~ 25), as done in the one- 
dimensional analysis by B05a. We also note that the dis- 
persion for g is modified as (j[g{0)] k, ag{\ — \g{p)\'^) in 
the subcritical, non- linear regime (Schneider et al. 2000); 
however, we neglect this non-linear correction for the dis- 
persion (jg to simplify various calculations, as in B05a. 
We found these are indeed good approximations for our 
combined distortion and depletion datasets, and the re- 
sults presented here are little changed when the full like- 
lihood function is used. 

The maximum likelihood solution, p, is obtained by 
minimizing the function F{p) with respect to p for given 
a and m. To do this, we compute numerically the 
derivatives dF{p) / dpi using a conjugate-gradient algo- 
rithm (Press et al. 1992). We determine by iteration the 
Bayesian value of a = d(p, m) by the following equation 
(Bridle et al. 1998): 



-2dS'(p, m) 



dTr(M-i) 



oodj 



(31) 



where M is a Afpix x A^pix matrix defined by 

M = G-i/=^HG-i/2^ (32) 



with 



d^F{p) 



d^Sip) 



dpidpj 



(33) 
(34) 



evaluated at p = p; iVgood is a measure of the effec- 
tive number of parameters (Suyu et al. 2006). With 
the optimal a, we thus expect that the final value of 
the misfit statistic is close to the classical number 
of degrees of freedom (hereafter NDF, see Suyu et al. 
2006), NDF = iVdaia - A^good (classic MEM).' We take 
TOj = const = m, where m can be regarded as a charac- 
teristic amplitude of the image (Maisinger et al. 1997). 
We find that the maximum-likelihood solution p for the 
Bayesian a is insensitive to the choice of m (see §5.5). In 
the following we set m to be 0.5. In order to be able to 



quantify the errors on the mass reconstruction we evalu- 
ate the Hessian matrix H of the function F{p) at p = p, 
from which the covariance matrix of the reconstructed k 
map is given by 

(35) 



Cij 



{6Ki5Ki) = {H)T.\p). 



4.4. Strong Lensing Constraints 

The map-making method described in §4.3 is only ap- 
plicable to subcritical regions lying outside of the critical 
curves (see Bartelmann & Schneider 2001). Therefore, 
our Subaru distortion and depletion datasets cannot be 
used to constrain the cluster mass distribution in the 
strong lensing region within the tangential critical curve 
(see B05b). Instead, we can pause quadratic constraints 
on such pixels to the total log-likelihood function based 
on strong lensing models as: 



F(p)-i^(p) + |^(«-«,)^ 



(36) 



where Kd is the constraint on Ki by strong lensing, Nc 
is the number of constraints, and /3 is a large number 
without causing this external term to become singular to 
working precision. 

In the present pixelization scheme, there is one such 
pixel (iVc = 1) containing the strong lensing region for 
which the Einstein radius of 9e — 43" corresponds to 
a mean redshift of — 0.87 for our red background 
population (see §3.3). To constrain the central k pixel, 
we utilized a mass model of A1689, which is well con- 
strained by ACS strong lensing observations restricted 
to the central region < 2' (B05b). The central ACS- 
derived mass profile of B05b is best described by an 
NFW model with a virial mass Mvir = 2.6 x IO^^Mq/H 
and a concentration Cyir = 8.2 having a scale radius 
of r, = rvir/cvir = 310kpc//i (B05b). With this NFW 
model, we find Kc ~ 0.781 ±0.1 at the central pixel of 
6 ~ (0', 0') for a reference source redshift of Zg = 1, cor- 
responding to Kc ~ 0.700 for our red background sam- 
ple (§3.3). Hence, assuming — 1 for the source red- 
shift will overpredict the central density by a factor of 
0.772/0.693 ~ 1.113, which will affect slightly the overall 
amplitude of the reconstructed k map (see §5.5). 

4.5. Downweighting Intrinsic Clustering Contributions 
and Noisy Measurements 

In contrast to the shearing effect, the magnification 
bias is a local, direct measure of the lensing convergence, 
free from the mass-sheet degeneracy. However, a prac- 
tical difficulty of the magnification bias measurement is 
the intrinsic clustering contribution which locally can be 
larger than the lensing induced signal in a given pixel. In 
order to obtain a clean measure of the lensing signal, such 
intrinsic clustering needs to be eliminated. Broadhurst 
et al. (1995) proposed an active declustering method 
based on the facts that the magnification bias depends 
strongly on the shape of the luminosity function, whereas 
intrinsic clustering depends weakly on the intrinsic lu- 
minosity function, and that they have different redshift 
dependence (see also Zhang & Pen 2005). This method 
however requires the addition of redshift information for 
individual background galaxies. Furthermore, the shape 
information of luminosity functions for respective back- 
ground populations must be provided in order to convert 
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a density depletion or enhancement into the lens magni- 
fication, /i. 

Alternatively, one may employ prior information from, 
for example, the surface luminosity density of cluster 
member galaxies to predict the projected mass distribu- 
tion, which is then used to downweight intrinsic cluster- 
ing contributions. In the present study, we have adopted 
an objective rejection scheme based on the gravitational 
shear predictions, summarized in the following steps: (1) 
Using shape distortion data alone we derive an entropy- 
regularized maximum likelihood solution p for the mass 
distribution. (2) Then, the shear-based k map is used 
to predict the magnification bias on the same grid of 
pixels. (3) Finally, we make a pixel-to-pixel comparison 
between the observed and predicted galaxy counts, and 
reject those magnification measurements Ni which are in 
conflict with the shear-based predictions Ni{p) as: 

\Ni - Ni{p)\ > ^/clust^/A^ (37) 

where A'o ~ 25 is the unlenscd, mean number counts of 
red galaxies, and fciust is a rejection threshold in units of 
a. We set a rejection threshold at 4(t, i.e., I'ciust = 4. 
Changing this threshold will affect the details of the 
reconstruction especially in the outer regions {9 ^ 7'). 
However, we confirmed that changing the threshold I'ciust 
by ±1 leaves our results unchanged within our statisti- 
cal uncertainties (see §5.5), and that the resulting mass 
map is fairly consistent with the surface luminosity and 
density distributions of cluster member galaxies, as de- 
scribed in detail below. 

In addition to the above, in order to exclude unreliable 
measurements, we have assigned zero weight to those pix- 
els which satisfy either of the following rejection criteria: 

1. Measurement pixel with no usable galaxy, 

2. Measurement pixel in the strong lensing region 

(§4.4), 

3. Magnification measurement near boundary regions 
with low completeness at fainter magnitudes. 

The last rejection criterion is required since the slope of 

the unlenscd number coimts, s = d\ogN()(rn)/dm, could 
depend on the selection completeness; such an apparent 
variation of the slope s could give rise to additional errors 
in the mass reconstruction. 

4.6. Map-making of A 1689 

We have applied our MEM method to the joint mea- 
surements of the shape distortion and magnification bias 
effects of red background galaxies in order to reconstruct 
the two-dimensional mass distribution of A1689. We uti- 
lized the FFTW implementation of Fast Fourier Trans- 
forms (FFTs) to convert between the lensing convergence 
and the gravitational shear fields using equations (7) and 
(9). The FFT however implies a periodicity in both hor- 
izontal and vertical directions, producing aliasing effects 
at the borders of the computational domain. In order 
to avoid such aliasing effects, we used large arrays of 
31 X 27 pixels with N^i^ = 31 x 27 = 837, covering a 
field of 43.'4 x 37.'8. On the other hand, the distortion 
and depletion measurements are limited to the central 
30' X 24' region (21 x 17 = 357 pixels); the wide-field 



Subaru data thus allow us to probe the projected mass 
distribution on scales ranging from 9 ^ 1' up to 6 ^ 20'. 
To minimize spurious effects from the periodic boundary 
condition, pixelized maps are further zero padded by a 
factor of 2 in each dimension (Seljak 1998; Sato et al. 
2003). 

For the spin-2 distortion measurements, we have in 
total A^data.ff = 2 X 355 = 710 usable (real) observa- 
tions. For the magnification measurements, we have 
.^data.ju = 302 usable observations. The total number 
of constraints is thus iVdata = A^daia.g + A^data.M = 1012, 
yielding A'data — Npi^ = 175 degrees of freedom (dof). 
The Baycsian value of the regularization parameter a 
that satisfies equation (31) was obtained as d = 96.6, 
when -2aS « A^pix - Q;Tr(M-i) « 213.1(= iVgood)- 
The minimum function value of F was found to be 
i^min = F{p) = {1/2)xHp) - aS{p) = 651.2, and 
X^(p) = 1089.3, corresponding to NDF = iVdata - 
Ng^od = 798.9, i.e., x^(p)/NDF = 1.36. Hence, the re- 
sulting is somewhat large, and is rather close to his- 
toric MEM, i.e., x^(p) ~ ^data- However, we note that 
our MEM mass reconstructions without the magnifica- 
tion data yield (p) /NDF 1 (see Table 2) as expected 
for classic MEM. Therefore, slightly large values of this 
misfit statistic using the magnification data could be at- 
tributed to the fact that the intrinsic clustering noise in 
the magnification measurements is not included in the 
likelihood calculation, underestimating the errors for the 
magnification measurements. 

In the left panel of Figure 7 we show the resulting k 
map on a grid of 31 x 27 pixels reconstructed from the 
combined distortion and depletion data with the ACS 
constraint on the mean value of k for the central pixel. 
The reconstructed spin-2 shear field (7[p]) is overlayed 
up on the k map. The right panel in Figure 7 shows 
the two-dimensional distribution of the rms reconstruc- 
tion error for Ki, estimated from the diagonal part of the 

pixel-pixel covariance matrix of errors: a{Ki) = C^/^{p). 
Figure 8 presents the contours of the reconstructed k map 
superposed on the V + i' pseudo-color image covering a 
field of 30' X 25' around the cluster. Here, for visualiza- 
tion purposes, the k map is resamplcd on to a finer grid 
and convolved with a Gaussian of FWHM = 1.'4, corre- 
sponding to a physical resolution of FWHM ~ 180kpc//i 
at the cluster redshift. The projected mass distribution 
of the cluster is smoothly varying and symmetric, with 
no significant substructure at r ^ 100kpc//i. Figure 9 
compares the reconstructed lensing fields and the mem- 
ber galaxy distributions in A1689, Gaussian smoothed to 
a resolution of FWHM = 2'. The top panels show the re- 
constructed K and the magnification-bias (5^ = ^^-^^^^^1 
fields in the left and right panels, respectively. The bot- 
tom left and right panels display the field-corrected i'- 
band luminosity and number density maps, respectively, 
of bright red-sequence galaxies with i' < 23 mag. For 
each panel the color scale is linear, and ranges from 0%- 
100% of the peak value. It is clear from Figure 9 that 
mass and light in A1689 are similarly distributed with 
a fairly round shape, and well centered on the main cD 
galaxy. 

4.7. Mass Maps from Different Datasets and Boundary 

Conditions 
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Any mass reconstruction technique based on the shear 
information involves a non local process, meaning that 
one has to assume certain boundary conditions to con- 
vert the gravitational shear field, 7(6), into the lensing 
convergence field, k{0) (Bartelmann & Schneider 2001; 
Umetsu et al. 1999). If the data field is sufficiently large 
so as to ensure that projected mass fluctiiations over the 
field average out, or that the mean convergence over the 
field is zero, then one may simply use equation (7) to 
invert the observed shear field into the convergence field 
(Kaiser & Squires 1993). Or, if magnification data are 
available, then a combination of complementary shearing 
and magnification information can be used to break the 
degeneracy between the obscrvables and the underlying 
gravitational potential field (Bartelmann 1995; Seitz et 
al. 1998). Besides, strong lensing observations, if avail- 
able, will place additional, tight constraints on the cen- 
tral mass distribution of the cluster where weak lensing 
alone cannot constrain (e.g., mass reconstruction of A370 
in Umetsu et al. 1999). 

Here we consider three sets of combinations of datasets 
and boundary conditions as summarized in Table 2: (i) 
2D MEM+ method using shear and magnification data 
with the central ACS constraint, (ii) 2D MEM method 
using shear and magnification data without the ACS con- 
straint, (iii) 2D MEM-S method using shear data with 
the central ACS constraint. For each MEM reconstruc- 
tion, we derive an entropy-regularized maximum likeli- 
hood solution for k with the Bayesian value of a. The 
resulting Bayesian value of a and the minimized func- 
tional values of F and are also listed in Table 2. 

In order to quantify the significance of the reconstruc- 
tion, we define an estimator for the detection signal-to- 
noise ratio (S/N) by the following equation: 

(S/N)2^ ^ ^ K^n.Cr-^, (38) 

«j>0 Kj>0 

where the indices i and j run over all pixels except those 
with negative values of n and those in the strong lens- 
ing region (see Table 2). The reconstructions based on 
both the distortion and depletion data yield a similar 
S/N of ~ 18-19, whereas the reconstruction from the dis- 
tortion data alone gives a slightly lower S/N of ^ 15 (sec 
Table 2). For comparison we quote the detection signifi- 
cance from the ID Subaru analysis by B05a based on the 
same red background catalogs for the weak lensing dis- 
tortion and depletion analyses: B05a measured the lens 
distortion and depiction profiles over a radial range of 
1' ;S 6* < 18'. We find S/N ~ 14.2 and S/N ~ 9.2 for the 
measurements of the distortion and depletion profiles, re- 
spectively. Since covariance between the distortion and 
magnification measurements can be neglected, the total 
S/N is simply obtained as \/14.22 + 9.2^ - 17. These 
numbers are quite comparable to those as measured from 
the reconstructed k field and its covariance matrix, indi- 
cating that the lensing information contained in the red 
catalogs is properly propagated into the Ac-basis'^. We 
find it important to use the Bayesian value for a in order 
to have an optimal smoothness for the mass reconstruc- 

This is not trivial since the noise level in a MEM-reconstructed 
map is affected by the smoothness constraint specified by a. Be- 
sides, MEM is non-linear, so that the resulting reconstruction er- 
rors depend on the signal as well as noise properties. 



tion, avoiding oversmoothing. A more quantitative com- 
parison between different reconstructions will be given in 
§5. 

5. MASS PROFILE OF A1689 

In this section, we aim to quantify and characterize the 
projected mass distribution (k) of A1689 reconstructed 
from Subaru weak lensing observations, in order to derive 
quantitative constraints on the three-dimensional mass 
distribution. Specifically, we will adopt the NFW den- 
sity profile p{r) cx r^^(l + r/rs)^"^ (Navarro et al. 1997) 
to describe the cluster mass distribution, characterized 
by the virial mass Myir and the concentration parameter 
Cvir = r^ii /rs, defined as the ratio of the virial radius Tvir 
to the scale radius rg . The best-fitting NFW parameters 
(Mvir,Cvir) can be then compared with ACDM predic- 
tions based on iV-body simulations (Bullock et al. 2001; 
Neto et al. 2007). As a test for the consistency and relia- 
bility, we will compare the best-fitting NFW parameters 
derived from different combinations of datasets, bound- 
ary conditions, and weak lensing techniques. Here we in- 
troduce three different methods to derive a convergence 
profile k{6) from weak lensing observations. 

5.1. Method (I): 2D Convergence Map 

The first method makes a direct use of the 2D k map 
reconstructed by the entropy-regularized maximum like- 
lihood method (§4). The k map is directly compared 
with the model convergence field for an NFW spheri- 
cal halo specified by two model parameters. We take the 
virial mass Mvir and the concentration parameter Cvir for 
describing an NFW halo. We employ the radial depen- 
dence of the convergence profile for the NFW model given 
by Bartelmann (1996). Note that the Bartelmann's for- 
mulae for the NFW convergence and shear profiles are 
obtained assuming the projection integral to extend to 
infinity. Alternatively, a truncated NFW profile can be 
used to model the convergence profile (Takada & Jain 
2003). We have confirmed that the best-fitting NFW 
parameters obtained using the above two different mod- 
els agree to within 1% for the case of A1689 lensing (for 
detailed discussions, see Baltz, Marshall, & Oguri 2007 
and Hennawi et al. 2007). We thus simply use Bartle- 
mann's formulae to calculate the relevant lensing fields 
for an NFW halo as done in B05a and B05b. 

We constrain the two NFW parameters from fit- 
ting to the 2D convergence map k{9) derived from Sub- 
aru weak lensing observations. We adopt a flat prior of 
Cvir < 30 for the halo concentration because the NFW 
profiles with Cvir ^ 20 cannot be distinguished by the 
Subaru data alone due to lack of information on the inner 
density profile (B05a). The x^-function for the Subaru 
weak lensing observations, x^ = Xwlj can be expressed 
as (Oguri et al. 2005) 

XwL = E [^(^') - '^(^')] - '^(^^•)] ' 

(39) 

where k{di) is the model prediction of the NFW halo for 
the ith bin {i = 1,2, ...,Nsl), and (C~^)y is the inverse 
of the pixel-pixel covariance matrix; TVpix = 31 x 27 = 
837. In the model fitting we exclude the central pixel in 
the strong lensing region (see §4.4). 
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We have derived sets of best-fitting NFW parameters 
for the three different MEM reconstructions listed in Ta- 
ble 2. Table 4 shows a summary of the best- fitting NFW 
parameters (Mvir, fVir) and the resulting value for our 
Subaru weak lensing observations; the errors quote 68% 
confidence intervals estimated from Ax^ = X^ — Xmin = 1 
in the (cyir, Mvir) -plane. Combining the distortion and 
magnification measurements with the central ACS con- 
straint (MEM-I-) yields the best-fitting NFW parameters, 
Mvir = (1.97 ± 0.20) X 10^^ Mq and Cvir = l^Atft, with 
Xmin/dof = 332/834(421), where the value in parenthe- 
ses refers to an effective degrees of freedom excluding 
upper limit bins with k < 0. As a test for the con- 
sistency, we compare best-fit NFW parameters for the 
three MEM reconstructions described in §4.7. Firstly, 
all of the three MEM reconstructions yield consistent 
results on the concentration parameter in the range, 
Cvir = 12.6 — 13.6, but with rather large uncertainties al- 
lowing a wide range of the concentration (8 ^ Cvir ^ 20). 
On the other hand, the best-fit values for A'/vir range from 
1.47 X lO^'^M© for MEM-S, through 1.62 x W^Mq for 
MEM, to 1.97 X IO^^Mq for MEM+, while the la error 
level for each is o-(Mvir) ~ 0.2 x IO^^Mq. 

The observed location of the Einstein radius can be 
used for a powerful, model independent test of the k pro- 
file (B05b). This is based on the fact that the enclosed 
mass interior to the Einstein radius is given by fundamen- 
tal constants and with knowledge of distances involved, 
namely k{9e) = 1 or M(< 6'e) = TT{Dd0E)'^^crit, pro- 
vided that the critical curve is nearly circular; this is the 
case for A1689 (B05b). Although our weak lensing mea- 
surements do not resolve such strong lensing phenomena, 
the observed Einstein radius of 9e = 45" (for Zs = 1) 
may be used to test the derived NFW models. To do 
this we numerically solve the equation 1 = knfw(^e) for 
the Einstein radius 9e by the Newton-Raphson method. 
We found 0E = 36.7;i8, 36.7^20, for the MEM-S, 
MEM, and MEM-I- reconstructions, respectively (Table 
4). All of the MEM reconstructions are roughly consis- 
tent with the ACS measurement of the Einstein radius; 
however, the constraints on the 6'e placed by the Subaru 
data alone are still rather weak, naturally due to the lack 
of central mass distribution. 

To make a direct comparison with model predictions, 
we compute the convergence profile k{6) from a weighted 
radial projection of the 2D n map as: 

Krr^ = K{0„^)= ^ W.^KiO,) / ^ Wi„, (40) 
iGBinm iGBinm 

where Wim (0 < Wim < 1) is the fraction of the area of 
the ith pixel lying within the mth annular bin. We use 
Monte Carlo integration to calculate these area fractions 
for individual pixels (Marshall et al. 2002), and 9m is 
the area-weighted center of the mth radial bin: 

9m= w^m/ j2 (41) 

i(EBinm iGBinm 

Since the Km profile is expressed as a linear combination 
of K{6i) values, it is straightforward to calculate the co- 
variance matrix Cmn = {5>^m^i^n) of the reconstruction 
errors: 



+ Y. W^mK-«(0O]y E (42) 

iGBinm \zGBinm / 

where the first term represents the errors propagated 
from the pixel-pixel covariance matrix of the 2D n re- 
construction, and the second term is responsible for vari- 
ations of the K(0)-field along the azimuthal direction, i.e., 
departure from circular symmetry. 

Here we derive for each of the 2D k maps a discrete 
convergence profile (40) in A^bin logarithmically spaced 
bins for 9 = [0min,0max]; wc adopt the same binning as 
in B05a: A^bin = 10, ^min = 1', and ^max = 18', as sum- 
marized in Table 3. In Figure 10 we compare the con- 
vergence profiles from the 2D MEM+ and MEM recon- 
structions, both of which are based on the lens distortion 
and magnification measurements, to clarify the effect of 
the ACS constraint on the central surface mass density. 
The error bars represent the \a uncertainties from the 
diagonal part (Cmm) of the bin-to-bin covariance ma- 
trix given by equation (42), and hence are correlated be- 
tween the different bins. One can see that overall the 
convergence profile obtained with the central ACS con- 
straint is slightly steeper than that without the ACS con- 
straint, and has a slightly higher overall normalization; 
these fcatiires are well explained by slightly higher NFW 
mass (Mvir) and concentration (cvir) derived for the re- 
sults with the central ACS constraint. Nonetheless, the 
two convergence profiles arc overall in good agreement 
within the statistical uncertainties. Steep NFW profiles 
with Cvir ^ 13 and Afvir ~ 2 x lO^'''A'/0 are well fitted to 
the reconstructed convergence profiles (solid and dashed 
curves). Figure 11 shows convergence profiles from the 
three MEM reconstructions: MEM-|- (squares), MEM 
(triangles), and MEM-S (crosses). Note that the vertical 
axis is linear here rather than logarithmic. The conver- 
gence profile from the distortion data alone (MEM-S) 
shows a slight negative dip of k; —0.01 at 9' ^ 6' due 
to boundary effects. 

5.2. Method (II): ID Maximum Likelihood Analysis 

In B05a we developed a method for reconstructing the 
discrete convergence profile Hm from a maximum like- 
lihood combination of radial profiles of the lens distor- 
tion and magnification effects on red background galaxies 
(see Figure 3 of B05a). With the assumption of quasi- 
circular symmetry in the projected mass distribution (or 
the projected potential field), one can express the lensing 
observables (i.e., tangential distortion and magnification 
bias) in terms of the binned convergence profile (see also 
§5.3). Owing to the nature of shear-based reconstruction, 
boundary conditions need to be specified (as discussed in 
§4.7). In B05a we set the inner boundary condition on 
the mass interior to ^min = 1', which is readily obtained 
from the well-constrained ACS mass model. 

One of the advantages of such a ID method is that 
one can improve up on the statistical significance of 
each measiircmcnt pixel by azimuthal averaging, pro- 
vided that the system is nearly symmetric. B05a mea- 
sured from the red background sample a tangential dis- 
tortion and a radial depletion profile in 10 discrete ra- 
dial bins, with total significance of S/N ~ 14.2 and 
S/N = 9.2, respectively (see §4.7). Thus, the per-pixel 
S/N of each measurement is of the order of unity, which 
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is optimal for an inversion problem. Accordingly, the ID 
analysis does not require regularization techniques. The 
best-fitting NFW parameters for the Subaru ID analysis 
are listed in Table 2. 

5.3. Method (III): Aperture Densitometry 

The spin-2 shape distortion of an object due to gravita- 
tional lensing is described by the complex reduced shear, 
g — gi + ig2 (sec equation [10]), which is coordinate de- 
pendent. For a given reference point on the sky, one 
can instead form coordinate-independent quantities, the 
tangential distortion g+ and the 45° rotated component, 
from linear combinations of the distortion coefficients gi 
and g2 as 

g+ = - {gi cos 2(l)+g2 sin 2^) , g^ = - (52 cos 2(j)-gi sin 20) , 

(43) 

where cj) is the position angle of an object with respect 
to the reference position, and the uncertainty in the g+ 
and fifx measurement is 0-+ = CTx = crg/-\/2 = cr in terms 
of the rms error (jg for the complex shear measurement. 
In practice, the reference point is taken to be the cluster 
center, which is well determined from symmetry of the 
strong lensing pattern. To improve the statistical signif- 
icance of the distortion measurement, we calculate the 
weighted average of the 5+'s and its weighted error as 



(^)=(Ev-0" 



■1/2 



(44) 
(45) 



For a shear-based estimation of the cluster mass profile 
we use a variant of the weak lensing aperture densitom- 
etry, or so-called the ^-statistic (Fahlman et al. 1994; 
Clowe et al. 2000), of the form: 



+- 



1 — (ftnn/^out)^ 
= k{6) - K(6li„„ < 1? < 6»out), 



dine' ^+{6') 



(46) 



where k{6) is the azimuthal average of the convergence 
field K{d) at radius 9, k{9) is the average convergence 
interior to radius 9, ^inn and ^out are the inner and outer 
radii of the annular background region in which the mean 
background contribution, K{, = R{9ian < < ^out), is de- 
fined; the "f+{9) = k{9) — k(9) is an azimuthal average 
of the tangential component of the gravitational shear 
at radius 9 (Fahlman et al. 1994), which is observable 
in the weak lensing limit: 7+(0) ~ {g+{9)). This cu- 
mulative mass estimator subtracts from the mean con- 
vergence k{9) a constant Kbg for all apertures 9 in the 
measurements, thus removing any DC component in the 
control region 9 = [^inn,^out]- We note that the Rb is a 
non-observable free parameter, and we use this degree-of- 
freedom to fix the outer boundary condition, and hence 
to derive a convergence profile. 

We compute the aperture densitometry profile Cc(^rn) 
in iV^ = 10 logarithmically spaced bins for 9 = 

[6*111111, ^'max], which we set to 6lniin = 1' and 6'max = 

16'; the maximum radius 6'max is comparable to the 
angular size of the cluster virial radius, ^vir = 15^7 



(rvir ^ 2Mpc//t), according to the ACS 4- Subaru- ID 
model (B05a). The inner and outer background radii arc 
set to ^inn = ^max = 16' and ^out = 19', respectively. Us- 
ing the ACS-|-Subaru-lD model by B05a, the mean back- 
ground level, Kb, is calculated to be Rb ^ 4.0 x 10~^ for 
our red background population with (Dds/Ds) = 0.693, 
at an effective source redshift oIzs^d = 0.68 (see equation 
[18]). Note that the current la upper limit is k{9) ^ 0.01 
at ^ ^ 8' (B05a). For a given boundary condition Rb, the 
average convergence R{9m) is estimated as 

R{0m) = Cc{Om) + Rb- (47) 

Then, we define a discretized estimator for k as 

Km = K{9m) = <C(^m+l) ' <C(^m) + « ' Oc'^Y^b, 

(48) 

where 



1 



2Aln6»™ V 6', 



1 



2Aln6'„ 



7m+l 



(49) 

and 9m is the weighted center of the mth radial bin 
(m = 1,2, ...,A^^ — 1; sec Appendix A). It is worth 
noting that, unlike strong-lensing based boundary condi- 
tions (e.g., B05a), this method utilizes an outer boundary 
condition on the mean background density Rb to derive 
a K profile (see Schneider & Seitz 1995 for an alternative 
method for a direct inversion of the mass profile). The 
error covariance matrix Cmn of Km is expressed as 



Crr 



-"l Q!2Wn,„+l - "2 "l <-m+l,„(50) 



C 



where = {6C,m5C,n) is the bin-to-bin error covari- 

ance matrix of the aperture densitometry measurements 
which is calculated by propagating the rms errors cr+(^m) 
for the tangential shear measurement (Okabe & Umetsu 
2008). 

In the non- linear regime, however, the 7+ {9) is not a di- 
rect observable. Therefore, non-linear corrections need to 
be taken into account in the mass reconstruction process. 
In the subcritical regime (i.e., outside the critical curves), 
the 7+(0) can be expressed in terms of the the averaged 
tangential reduced shear as {g+{9)) w ^+{9)/[1 — k(9)] as- 
suming a quasi circular symmetry in the projected mass 
distribution (B05a; Umetsu et al. 2007). This non-linear 
equation (46) for Cc(^) can be solved by an iterative pro- 
cedure: Since the weak lensing limit (k, I7I, \g\ <C 1) holds 
in the background region 6'inn ^ 9 < 6'niax, we have the 
following iterative equation for Cc(^): 



(fc+i; 



{9)^2 j '"°d\n9'{g+{9))[l-K^ 
Je 

ut 

d\n9'{g+{9')),{bl) 



1 — {9inn/9ontY 



where Cc*^^^' represents the aperture densitometry in the 
{k + l)th step of the iteration (fc = 0, 1, 2, Mter); the 
^(fc-i-i) Jg calculated from ^c'^'''^^ using equation (48). This 
iteration is preformed by starting with k*^^) = for all 
radial bins, and repeated until convergence is reached at 
all radial bins. For a fractional tolerance of 1 x 10~^, 
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this iteration procedure converges within iViter ~ 10 it- 
erations. We compute errors for (c and k with the hnear 
approximation. 

Figure 12 shows the resulting model-independent mass 
profile (squares) of A1689 with decorrelated error bars, 
along with the results without the non-linear corrections 
(triangles). Without the non- linear corrections, central 
bins at 9^3' are underestimated by 15% at maxi- 
mum. The reconstructed convergence in the first bin is 
k{1'.2) = 0.32 ± 0.15 (decorrelated Icr error), which is 
consistent at the ~ la level with the previous Subaru 
ID results (B05a): ^(1^2) = 0.44^°;^^ (Icr) for the red 
background population with {Dds/Ds) ~ 0.693. We note 
that B05a utilized an inner boundary condition on the 
mean k interior to ^min = 1' based on the ACS mass 
profile. 

Force fitting an NFW profile to the derived k profile 
yields Mvir = (1.48±0.27) x IO^^Mq and Cvir = '^T-^tHa 
with the adopted prior Cvir < 30, where x^jj^/dof = 
5.2/7(7). It is interesting to compare the reconstructed 
K profile with the tangential shear profile, (g_|_(^)), which 
is a direct observable with uncorrelated measurement er- 
rors (see equation [44]). In B05a we show the tangential 
distortion profile for the same red background sample as 
used in the present study (see Figure 1 of B05a, or Fig- 
ure 16 below). The best-fitting NFW model for the shear 
profile is given by Mvir = (1.51 ± 0.26) x IO^^Mq and a 
high concentration, c^i, = 20tli (Xmin/dof = 5.0/8), 
which is in good agreement with the results from the 
ID reconstruction based on the aperture densitometry. 
Such a detailed agreement ensures the validity of the 
boundary condition for a shear-based mass reconstruc- 
tion. We note that simply assuming Kb = yields sim- 
ilar results, Mvir = (1-51 ± 0.3) x IO^^Mq (cvir < 30) 
with Xmin/'lof = 5.5/7(5), being consistent within the 
Icr uncertainty. 

5.4. Combining Strong and Weak Lensing Results 

ACS Constraints 

As demonstrated in B05a and the previous sections, it 
is crucial to have information on the central mass distri- 
bution in order to derive useful constraints on the con- 
centration parameter, and hence to distinguish the NFW 
profile from others. To do this, we constrain the two 
NFW parameters from fitting to the combined Sub- 
aru weak lensing and ACS strong lensing observations 
following B05a and Oguri et al. (2005): 

X'=XWL+XSL, (52) 

where the Xwl ^'^^ the Subaru observations is defined 
by equation (39). For the ACS data we use an az- 
imuthally averaged profile of the k map well constrained 
by the strong lensing observations of B05b; this profile 
is given in A^sl = 12 bins linearly spacing over the range 
e = [0^077, 0f97] (see Figure 22 of B05b, or Figure 3 of 
B05a), and the amplitude is scaled to (Dds/Dg) = 0.693 
of our red backgroimd sample from Dds/ Dg = 0.881 at 
Zs = 3. It is important to note that the ACS strong lens- 
ing analysis of B05b unveiled the secondary mass clump 
associated with a small clump of galaxies (Teague et al. 
1990; Czoske 2004). However, the mass contribution of 
this subclump is only a small fraction of the cluster mass 
component (see Figures 21 and 22 of B05b) and it only 



slightly perturbs the tangential critical curve (see B05b). 
The ACS mass profile is corrected for the subclump, 
where the secondary clump region is locally masked out 
when taking an azimuthal average (B05b). By combining 
the Subaru and ACS lensing analyses, we can trace the 
cluster mass distribution over a large range in amplitude 
K ^ [lO^'^, 1] and in projected radius r = [10^^, 2]Mpc//i. 
The function for the ACS observations is expressed as 

= E ^^^^2^ (53) 

where k{Gi) is the model prediction of the NFW halo for 
the ith pixel {i = 1,2,..., A/pix), and ai is the Icr error for 
Kf. the bin width of the ACS-dcrived convergence profile 
is sufficiently broad to ensure that the errors between 
different bins are independent (B05b) 

Einstein- Radius Constraint 

Alternatively, as a model independent constraint,* we 
can utilize the observed location of tangential critical 
curves traced by giant arcs and multiply imaged back- 
ground galaxies, defining an approximate Einstein ra- 
dius, ^E- This radius is determined from strong lensing 
modeling of many multiple images visible around these 
clusters in deep HST/ACS images (B05b), and corre- 
sponds to the theoretically extreme value of the elliptic- 
ity, g+{OE) = 1- This Einstein-radius constraint comple- 
ments in a model independent manner the weak lensing 
shape measurements. 

We constrain the NFW model parameters (cvir,Mvir) 
by combining weak lensing profiles and the Einstein- 
radius constraint. We define the x^ function for com- 
bined weak lensing distortion and Einstein-radius con- 
straints by 

2 _^ [{g+{em))-g+{0m)f . [l-g+{0B,ZE)f 

(54) 

where the first term is the x -function for the Subaru tan- 
gential shear measurements (§5.3) and the second term is 
the x^ function for the Einstein radius constraint; g+{Om) 
is the NFW model prediction for the reduced tangen- 
tial shear &i 6 = 9m calculated for the red background 
sample (see §3.3), g+(6'E, ^e) is the model prediction for 
the reduced tangential shear at ^ = ^e, evaluated at 
the arc redshift, Zg = Ze, and ct+^e is the rms error for 
g+{9E, Ze)- Following B05b, we take 9e = 45" at ze = 1, 
and assume conservatively a 10% error for the Einstein 
radius: (T0^/9e = 0.1. We then propagate this error to 
g+ assuming a singular isothermal sphere (SIS) model; 
At r ^Vs, the density slope of NFW is shallower than 
that of SIS (see Figure 1 of Wright & Brainerd 2000), 
so that this gives a fairly conservative estimate of cr_)_^E- 
For the SIS model, d\n g+/dhi9E = 2 at 9 = 9e, so that 
£7+,E = 2ae^/eE = 0.2. 

Alternatively we can combine the k profile recon- 
structed from weak lensing distortion measurements 
(§5.3) with the Einstein-radius constraint. For an axially 

* Wc note that this is based on the assumption of the circu- 
lar symmetry in the projected lens mass distribution, and that a 
rather tight prior (estimated from data) is applied to one of the 
"parameters" (i.e., 0^)- 
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symmetric lens, the Einstein-radius constraint is simply 
written as k(0e,^e) = 1- Thus, the function is now 
expressed as 



the non-linear Cc-statistic measurements. The filled cir- 
cles show the results from the Subaru ID analysis by 
B05a based on the combined magnification and distor- 
tion profiles of the same red background population as 
ii{Sm Hoigjpresent study. It is clearly seen from Figure 13 that 



(55) 

where k{0-e,ze) is the model prediction for the average 
convergence interior to the Einstein radius 9e evaluated 
at the arc redshift of Zg = Ze, and a^^E is the rms error 
for k{6e, ze), which we take as a^^E = c^e/^e = 0.1. 

Comparison 

Table 4 summarizes for each mass reconstruction the 
best-fitting NEW parameters (Mvir,Cvir), the minimized 
value with respect to the degrees of freedom, and 
the predicted Einstein radius 9e for a fiducial source at 
Zs = 1. Here for each of the Subaru reconstructions, we 
compare the best-fit NEW parameters with and without 
the ACS inner profile combined, as indicated in the third 
column of Table 4. When combined with the inner ACS 
profile, all of the Subaru reconstructions yield a virial 
mass in the range M^ir ^ (1.8 — 2.1) x lO^^M© and a 
high concentration in the range Cyir ^ 13 — 15, prop- 
erly reproducing the observed Einstein radius of ~ 45" 
at = 1 (B05b). In particular, fitting an NEW pro- 
file to the combined ACS and Subaru-2D data based 
on the MEM-I- mass reconstruction (see Table 2) yields 
Mvir = (2.10 ± 0.17) X IQi^Mo and = 12.7^0;" 
(Xmin/dof = 335/846; effective dof of 433 without in- 
cluding upper limit bins with «; < 0), with the predicted 
Einstein radius of 6e = 45.3^g'2 arcsec. 

The 2D-bascd results here can be compared with the 
corresponding results from the Subaru ID analysis of 
B05a (see Table 4) based on the combined distortion 
and magnification measurements. The combined ACS 
and Subaru-ID convergence profile is well fitted by an 
NEW profile with Mvir = (1.9 ± 0.2) x IO^^Mq and 
Cvir = 13.7lJ;^ (B05a), with the predicted Einstein ra- 
dius of Be = 45.4lg;g arcsec, in good agreement with 
the present full 2D analysis within the statistical uncer- 
tainties; this agreement between the ID and 2D anal- 
yses supports the assumption of quasi-circular symme- 
try in the projected mass distribution. It is interesting 
to compare these results with different combinations of 
lensing measurements having different systematics. An- 
other combination of the ACS and Subaru ID conver- 
gence profile, derived from the shear-based ^-statistic 
measurements (§5.4.0), yields fairly consistent results: 



M^w = 1.91 



+0.24 
-0.20 



x liV^Mc; and Cv 



13.7 



+ 1.5 
-1.3- 



This 



consistency clearly demonstrates that our results here 
are insensitive to systematic errors in the lensing mea- 
surements, such as the shear calibration error. 

In Eigure 13 we make a direct comparison between 
model-independent convergence profiles of A1689 from 
different weak lensing tcchniciues, along with the ACS- 
derived inner profile obtained by B05b (filled triangles). 
The filled circles with error bars show the results based on 
the 2D K map reconstructed from an entropy-regularized 
maximum-likelihood combination of the lens magnifica- 
tion and distortion of red background galaxies (MEM-I- , 
Table 2). The open triangles represent the k profile from 



(j|the Subciru-derived k profiles are all in good agreement 
within the statistical uncertainties, over the full range of 
radii up to ~ 2Mpc//i. The entire mass profile traced by 
the combined ACS and Subaru information is well de- 
scribed by a single NEW profile (solid and dashed) with 
a high concentration (cvir ~ 13— 14). Eor comparison an 
NEW profile with a low concentration Cvir ~ 5 normal- 
ized to the observed Einstein radius of 45" (at Zg = 1) 
is shown as a dotted curve. Such a low concentration 
profile as favored by the standard ACDM model clearly 
overpredicts the outer profile constrained by the Subaru 
weak lensing observations. 

Eigure 14 shows the 68%, 95%, and 99.7% confidence 
levels (Ax^ = 2.3, 6.17, and 11.8) in the (cvir, Mvir)-plane 
for each of the three-different 2D MEM reconstructions 
(from top to bottom: MEM-S, MEM, and MEM-F, see 
Table 2) with (right) and without (left) the central ACS 
profile at 10kpc//i ^ r ^ 180kpc//i combined. As shown 
in the left panels, Mvir is well constrained by the Sub- 
aru data alone, while the Subaru constraint on Cvir is 
rather weak, allowing a wide range of the concentra- 
tion parameter, Cvir. The complementary ACS obser- 
vations, when combined with the Subaru observations, 
significantly narrow down the uncertainties on Cvir (right 
panels), placing stringent constraints on the inner mass 
profile. For each sub-panel, the observed constraints on 
the Einstein radius, — 45" at z,; = 1, are shown as a 
dashed curve. This clearly demonstrates that the ACS- 
derived constraints on the K-field ensure the correct size 
of the observed Einstein radius. 

Table 5 lists the best-fitting NEW parameters obtained 
in different studies. Here wc quote as the NEW param- 
eters (c2oo,.^20o) evaluated at a specific fractional over- 
density of Ac = 200 with respect to the critical density 
Pciit{zd) for closure of the universe at Zd = 0.183 as well 
as those in terms of the virial properties, (Cvir, Mvir). We 
convert a given set of (c2oo, M200) into (cvir, Mvir) assum- 
ing a spherical NEW density profile (see, e.g.. Appendix 
A of Shimizu et al. 2003). We show in Figure 15 the same 
confidence levels as in Eigure 14 but for the Subaru g+ 
profile of B05a. Also shown are best-fit sets of (cvir, Mvir) 
taken from Halkola et al. (2006) and Limousin et al. 
(2007) as well as from our combined ACS and MEM-I- 
results. The best-fitting NEW parameters for the Sub- 
aru g+ profile (cross) are Mvir = 1-51^0 24 ^ 10"'^''' M0 and 
Cvir — (Table 4), being marginally consistent 

with the very high concentration of Cvir = 27.2^5;^ de- 
rived by Medezinski et al. (2007) ^. The k field derived 
from the present full 2D analysis, in conjunction with the 
central ACS profile, favors a slightly lower, but still high, 
concentration of 13 (triangle), reproducing the ob- 

served Einstein radius of 45" at Zg = 1. On the other 
hand, Limousin et al. (2007) found from their CFHT 
weak lensing analysis a concentration of Cvir = 9.6 ± 2.0 

^ Unlike our analysis here, Medezinski et al. (2007) used the 
observed constraints on the Einstein radius when fitting an NFW 
profile to the g+ profile measured from the combined red and blue 
sample of the background. 
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and a virial mass of Mvir = l-StoJ ^ lO^^Mg (filled 
circle), corresponding to the Einwtcin radius of 24 ± 11 
arcsec (at Zs = 1), inconsistent with the observed Ein- 
stein radius. The results from Halkola et al. (2006) 
with Cvir = 9.6 ± 0.4 and M™ = ^.Gtal x W^^Mq 
(square) come closer to the observed Einstein radius, 
being marginally consistent within the la uncertainty 
{6e = 39^3 arcsec at Zg = 1), based on almost exactly 
the same shear data used in B05a but with a different 
weighting which prefers their inner strong-lensing based 
profile where the ACS data imply a shallower slope, as 
discussed in Medezinski et al. (2007) and Limousin et al. 
(2007). 

It is useful to compare the results from different lens- 
ing studies in terms of the tangential distortion, that 
is directly observable in weak lensing. Figure 16 shows 
the radial profiles of g+ and gx derived in B05a from the 
same red background sample as used in this work, where 
we have added the observed Einstein-radius constraint 
of (?E = 45 ± 5 arcsec {zg ~ 1), translated to the mean 
depth of the red backgroimd sample (§3.3), marking the 
point of maximum distortion, = 1 {0^ = 39 ± 4 arc- 
sec). The ACS-|-Subaru-2D NEW model (sohd curve) 
fits reasonably well with the entire distortion profile, 
r ^ [80, 2000]kpc/ft,, although it slightly overpredicts the 
outer profile at 9 ^9', meaning that the observed g+ pro- 
file is steeper at large radii, as pointed out in Medezinski 
et al. (2007). In contrast the best-fit NEW model of 
Limousin et al. (2007), shown with the dashed curve, is 
in agreement with the Subaru outer profile, particularly 
at 10' ^ 9 ^ 18', but underpredicts significantly the in- 
ner g^ profile, leading to a significant underprediction 
of the Einstein radius. On the other hand, the dotted 
curve in Figure 16 shows the NEW profile of Halkola et 
al. (2006), which is overall in good agreement with the 
observed g+ profile, and with the NEW prediction of this 
work, but increasingly overestimates the distortion signal 
with radius, and slightly underpredicts the Einstein ra- 
dius. This is again consistent with that the discrepancy 
between the derived NEW parameters of this work and 
Halkola et al. (2006) is due to the relative weights in the 
fitting procedure assigned differently to the strong and 
weak lensing measurements; that is, a relatively higher 
weight is given by Halkola et al. (2006) to the shallower 
inner profile at r^40kpc//i constrained by the radial 
arcs,^° and consequently, the location of the outer criti- 
cal curve is slightly underpredicted. However, we argue 
this does not indicate a discrepancy between our strong 
and weak lensing results, but is simply that the form of 
the NEW profile is not entirely consistent with our data 
ranging from 10 to 2000kpc//i; the best-fit NEW profile 
is either too shallow at large radii or too steep at small 
radii, depending on the radial limits being examined. 

5.5. Systematic Errors of the Concentration Parameter 

In this subsection we address the issues of systematic 
errors on the halo concentration Cvir inherent in our lens- 
ing measurements and analysis methods. Here we dis- 
cuss the following potential sources of systematic uncer- 
tainty: (i) selection criteria for the red background sam- 
ple (§3.3), (ii) background redshift distribution (§3.3), 

The location of the observed radial critical curve is at ~ 17" 
as found in B05b. 



(iii) magnification analysis (§4.2, §4.5), (iv) inner bound- 
ary condition (§4.4), (v) strong lensing model (§4.4 and 
§5.4.0), (vi) entropic prior (§4.3), and (vii) shear cali- 
bration error. In Table 6 we summarize the systematic 
errors of the halo concentration, Cvir, given in fraction of 
the best-fit value of Cvir = 12.7 obtained from the com- 
bined ACS and Subaru-2D (MEM+) data. Adding all 
potential sources of error in quadrature, the total uncer- 
tainty in our determination of the halo concentration is 
Cvir = 12.7 ± 1 (statistical) ± 2.8(systematic). 

Background Selection Criteria 

As clearly demonstrated by B05a and Medezinski et 
al. (2007), a secure background selection is critical in 
the cluster weak lensing analysis, in order to avoid dilu- 
tion of the distortion signal by contamination of unlensed 
cluster member galaxies. The degree of dilution, which 
is proportional to the fraction of cluster membership, is 
particularly prominent in the central region of rich clus- 
ters, as in the case of A1689 (see Medezinski et al. 2007). 
Practically, a reliable background sample can be defined 
by selecting galaxies with colors redder than the cluster 
E/SO sequence (see §3.3). Eigure 4 shows that the di- 
lution effect becomes significant when the lower (bluer) 
color limit of the entire red sample is decreased below a 
color of A{V — i') ~ 0.1, while no signature of systematic 
variations is seen when the lower color limit is increased 
above A{V — i') ^ 0.1, ensuring that the dilution effect 
is almost negligible there. Based on this, we defined our 
red background sample by choosing a conservative color 
limit as A(F - i') > 0.22. 

Here we vary the lower color limit of the red sample in 
the interval of [0.1,0.5] where the dilution effect is neg- 
ligible. At the lower color limit of ^ 0.5, however, the 
number of red galaxies is decreased by about 40% (see 
Eigure 4), when compared with our fiducial sample of 
A{V — i') > 0.22, leading to noisier results. We gen- 
erate red samples at the lower color limit of 0.1, 0.15, 
0.2, 0.5, and assess the scatter in the best-fit concen- 
tration parameter, Cvir- We find the error distribution 
is fairly random, with a small spread of Acvir ~ ±0.5, 
corresponding to a ~ 4% fractional systematic error. 

Background Redshift Distribution 

Redshift information of background galaxies plays a 
crucial role in the determination of cluster mass profiles 
In order to convert the observed lensing signal into phys- 
ical mass units, one needs to evaluate the mean distance 
ratio (Dds/Dg) over the redshift distribution of back- 
ground galaxies (§3.3). An overestimate of the source 
redshift will systematically lead to an underestimate of 
the cluster mass, and vice versa. In this way, the uncer- 
tainty in the background redshift distribution will lead to 
systematic errors of the cluster mass determination. The 
level of uncertainty depends on the lensing geometry, and 
is less significant for low-z clusters (say, Zd ^ 0.2). Eor 
purely weak lensing based data, this effect is less impor- 
tant for the determination of the halo concentration, be- 
cause the identification of the inner characteristic radius 
(rg) is basically independent of the background redshift. 
However, when the weak lensing measurements are com- 
bined with inner strong lensing information, this depth 
information becomes crucial for the determination of halo 
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concentration as well, because it determines the relative 
normalization between the inner and outer profiles. That 
is, an overestimate of the background redshift will cause 
an underestimate of the surface mass density, which will 
increase the difference between the inner profile derived 
from strong lensing and the outer profile from weak lens- 
ing, leading to an overprediction of the halo concentra- 
tion. 

Here we turn to assess the level of systematic error 
arising from the uncertainty in the background redshift 
distribution. Based on the multicolor photometry of Ca- 
pak et al. (2004) in the HDF-N, we obtain a mean dis- 
tance ratio of {Dds/Ds) = 0.693 ± ±0.02, or a distance- 
equivalent redshift of Zg^jj = 0.68 ± 0.05, for onr color- 
magnitude selected red background sample (§3.3; see also 
Medezinski et al. 2007). A good agreement has been 
also found using the COSMOS deep multicolor photo- 
metric catalog of Capak et al. (2007), yielding a similar 
depth of (Dds/Ds) « 0.703 (Medezinski et al. 2007 in 
preparation), suggesting that the field- to-field variance 
(cosmic variance) is not significant, and as small as the 
statistical uncertainty obtained here. We find this level 
of depth uncertainty will lead to only a ~ 1% fractional 
systematic error of Cvir- Further, we assign a conservative 
uncertainty in the (distance-equivalent) source redshift, 
Zs,D = [0.7, 1.0], and find a fractional systematic error of 
about 10% in the concentration parameter. It is interest- 
ing to note that B05a obtained the best-fit concentration 
parameter of Cvir = 13.7lJ'i assuming a source redshift 
of Zs{= Zs,d) = 1- This concentration is slightly higher 
(~ 4-8%) than the best-fit value of Cvir = 12.7 ± 1 found 
in this work, but this level of discrepancy can be easily 
reconciled by the systematic error in the assumed depth 
of B05a: Applying this bias correction to the results of 
B05a yields Cyir ~ 12.3, which agrees quite well with the 
present results of Cvir = 12.7 ± 1. 

Lensing magnification influences the observed surface 
density of background galaxies (Broadhurst et al. 1995), 
by expanding the observed solid angle of the background 
(area distortion), and decreasing the effective fiux limit of 
the survey (flux amplification). Thus the latter effect of 
magnification bias may change the redshift distribution 
of background galaxies depending on the distance from 
the cluster center, which could be a potential source of 
the systematic errors in the determination of the halo 
concentration. Our red background sample, however, is 
highly depleted (see Figure 6), meaning that the area 
distortion effect is dominating over the fiux amplifica- 
tion of fainter, distant background galaxies. Indeed, the 
unlensed count slope measured at our magnitude limit 
i' = 25.5 is s = d log No {m) / dm « 0.22, being fairly flat 
as compared to the blue sample with s w 0.4, close to 
the lensing invariant slope (see equation [23]). Conse- 
quently, relatively few fainter objects are magnified into 
the sample even in the central cluster region, so that the 
magnification effect on the source redshift distribution is 
fairly negligible for the red background defined at fainter 
magnitude limits For our sample, magnification at 9 ^ 2' 
is ~ 1 -h 2k « 1.2, or about 0.2mag of increased depth, 
corresponding to the fractional increase of only ~ 5% in 
the number of red background galaxies. Further, given 
the weak dependence of the redshift distribution of faint 
galaxies on apparent magnitude (e.g., Medezinski et al. 



2007), we do not need to take this effect into serious con- 
sideration for our red background sample. 

Still, it is instructive to consider the potential system- 
atic bias caused by the magnification effect on the back- 
ground redshift distribution. The net effect of depth cor- 
rection is to reduce the central surface mass density of 
the lens, since we attribute the high lensing signal to 
the geometric information of the background. For purely 
weak-lensing based data, this will lead to a lower concen- 
tration. However, when the strong lensing information in 
the inner region is taken into account as well, this depth 
correction will further increase the difference between the 
inner and outer profiles derived from strong and weak 
lensing, respectively, thereby enhancing the concentra- 
tion. However, the amount of correction is negligible in 
practice when the weak lensing analysis is based on the 
red background galaxies as discussed above. 

Magnification Analysis 

Here we address the systematic uncertainties arising 
from a particular treatment and various cuts in the weak 
lensing magnification analysis. 

(1) Mask area correction 

The masking effect due to cluster member galaxies and 
bright foreground objects acts to reduce the apparent 
number of background galaxies, and this reduction in- 
creases towards the cluster center, leading to an overes- 
timate of the central depletion signal without the mask- 
ing correction. In the present study the mask area of 
bright objects is evaluated as the area inside the ellipse 
of fmask-times the major (A_IMAGE) and minor (B_IMAGE) 
axes computed from SExtractor (see Cobb et al. 2006 
for a similar discussion), where the multiplier is chosen 
as i^niask = 3 SO that the ellipse is visually consistent with 
the isophotal detection limit in our SExtractor configu- 
ration (§3.1). Here we adopt a conservative uncertainty 
of t'mask = 3 ± 1 on the masking factor, and find the cor- 
responding systematic uncertainty of ±7% in the halo 
concentration. The lower limit on Cvir is obtained for 
t'mask = 2, when the mask area correction is underesti- 
mated and hence the depletion signal is overestimated. 

(2) Clustering rejection 

Similarly, we adopt a conservative uncertainty of 
i^ciust = 4 ± 1 in the rejection threshold t'ciust) which 
has been introduced to downweight locally the intrinsic 
clustering noise which otherwise perturbs the depletion 
signal (§4.5). This yields a fractional uncertainty of only 
±1.5% in Cvir- We note that the intrinsic clustering of 
background galaxies is a local effect, and hence it does 
not affect significantly the radial profile fitting. 

(3) Unlensed count slope 

The conversion from the observed counts of the 

red background sample into the magnification bias 5^ 
(equation [23]) depends on the slope parameter s = 
d\ogNi)(rn) / dm, of the unlensed number counts (§4.2), 
which was estimated as s = 0.22 ± 0.03 from the outer 
region ^ 10' where the magnification effect is negligi- 
bly small (k, 5^ ;$ 0.01). We find that this level of un- 
certainty in s will cause an uncertainty in Cvir of about 
3.5%. 
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Inner Boundary Condition 

Any mass reconstruction technique based on the grav- 
itational shear field involves a non local process (see 
equation [6]), and hence it is crucial to have a proper 
boundary condition for an accurate determination of the 
cluster mass profile. In particular, our MEM+ method 
is based on the ACS strong lensing information, which 
is incorporated as an inner boundary condition on the 
central pixel, Kc (§4.4). B05a showed that the com- 
bined ACS and Subaru mass profile can be well fitted 
by a high concentration NFW profile (cyir ~ 14) over 
the full range of ACS and Subaru data, r = [10~^,2] 
Mpc/h. However, it is also found in B05a that this high- 
concentration model somewhat overestimates the inner 
slope at r ;S 40 kpc//i (sec Figure 13). This model yields 
a central surface mass density of Kc{zs = 1) « 0.810, 
which is slightly higher than, but still consistent with, 
the prediction based on B05b adopted in the present 
work, Kc{zs = 1) = 0.781 ± 0.1, within the la statisti- 
cal uncertainty. The mass profile of A1689 has also been 
examined by Limousin et al. (2007) using independent 
weak lensing shape measurements from CFHT12K data. 
Their best-fit NFW model, however, underpredicts the 
observed Einstein radius (§5.4.0), and accordingly yields 
a much lower central value of k,c{zs = 1) ~ 0.527, which 
we take as the lower limit on Kc{zs = 1). Allowing for a 
conservative uncertainty of Kc{zs = 1) = 0.78lo!25' 
responding to Kc(2s,r> = 0.68) = 0.70^0.22 for the effec- 
tive depth of our red sample, we find a ±10% fractional 
systematic uncertainty in the halo concentration, Cyir- 

Inner Strong Lensing Information 

The inner strong lensing information provides strong 

constraints on the halo concentration parameter as 
demonstrated in Figure 14. In the present work, the 
ACS-derived inner mass profile of B05b is specifically 
used to determine the NFW halo parameters of A1689 
in the strongly lensed region, r = [10,130] kpc//i. It 
is practically difficult to assess potential systematic er- 
rors introduced in strong lensing modeling because of 
the complex, non-linear error propagation (B05b). In- 
stead, here we simply estimate the level of uncertainty by 
a comparison with the result obtained with the model- 
independent constraint on the Einstein radius 9^ (see 
§5.4.0) based on multiply lensed images identified in the 
ACS observations of B05b. Both the strong lensing anal- 
ysis of B05b and Limousin et al. (2007) yield a con- 
sistent value for the projected mass interior to 45" of 
M2d(45") w 2 X IO^^^Mq, or equivalently k(45") = 1 for 
a source redshift of = 1: i.e., O-e = 45" at Zg = 1. In 
contrast to the fit to the inner mass profile, this Einstein- 
radius information provides an integrated constraint on 
the inner mass profile interior to 45", or r « 100kpc//i 
in projected radius. We find that a joint fit of the Sub- 
aru K map and the Einstein-radius constraint yields a 
slightly higher concentration of Cvir = 14.0t2'i (see Ta- 
ble 4), corresponding to a fractional increase of about 
10%. This tendency is also found for the results with the 
shear-based ID mass reconstruction from the Cc-statistic 
measurements (see Table 4). 

Entropic Prior 



A particular choice of the regularization could be a 
potential source of the systematic errors in the cluster 
mass reconstruction. In our mass reconstructions the 
model parameter m of the entropy prior is fixed, but the 
Bayesian value of the regularization parameter a that 
satisfies equation (31) is obtained for a given value of m. 
When we vary the value of m over a relevant range of the 
cluster lensing signal, m = [0.1,0.9] (m — 0.5 ± 0.4), we 
find the error distribution in the resulting value of CyJ J- IS 
almost random with a small spread of ±10%. 

Furthermore, a particular choice of the entropy func- 
tion may lead to some systematic bias in the determi- 
nation of cluster mass profiles. To check this possibil- 
ity, here we simply compare the present results from 
the MEM+ reconstruction with earlier ID maximum- 
likelihood results of B05a, both of which are based on the 
same distortion and magnification data, and adopt the 
ACS-based inner boundary conditions. After the correc- 
tion for the systematic bias (§5.5.0) the best-fit concen- 
tration of B05a is Cyir ~ 12.3, which is in good agreement 
with Cvir = 12.7 ± 1 obtained with the entropy regular- 
ization. Thus, it is likely that the level of systematic 
uncertainty due to the particular choice of the entropy 
regularization is negligibly small as compared to other 
sources of the systematic errors. 

Our use of the entropy prior in the low S/N regime 
might potentially induce some slight bias in the regu- 
larized maximum likelihood solution, and/or some slight 
non-Gaussianity in the error distribution, underestimat- 
ing the actual error bars for the mass reconstruction and 
the NFW halo parameters. However, our results show 
good consistency between the entropy-regularized recon- 
structions and other standard reconstructions within the 
statistical uncertainties. Thus, it seems this bias is not 
significant for this work. 

Shear Calibration Error 

A shear calibration error is one of the systematic errors 
that could bias the weak lensing shape measurements 
(Haymans et al. 2006; Massey et al. 2007) and poten- 
tially have some influence on recovered mass profiles. To 
assess this possibility we have measured the strength of 
the weak lensing signal as a function of magnitude limit 
for our background galaxy sample. We found no partic- 
ular tendency towards a loss of the weak lensing signal 
with apparent magnitude for red background galaxies, 
which span over a wide range of sizes. This is comfort- 
ing and consistent with the expectations of the model- 
independent KSB± based technique for which the recov- 
ered signal should match reality within the noise. We 
note that at the very weak lensing limit, our distortion 
measurements are quite consistent with an independent 
estimate of the weak lensing signal by Limousin et al. 
(2007). 

However we have found that for blue background galax- 
ies there is a significant loss of the signal at faint magni- 
tudes and this raises the worrying question of unresolv- 
ablc HII regions which we know do become prevalent at 
faint blue magnitudes, acting effectively as point sources 
and hence reducing the weak lensing signal. Such ob- 
jects are not included in our analysis, so as not to bias 
our lensing measurements. The STEP project, aimed 
at assessing signal which may be lost in ground based 
data, described in Heymans et al (2006) and Massey et 
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al. (2007), docs not allow for unresolvable sources within 
galaxies, as it is inherently assumed that galaxies are 
continuously resolvable, so that stretched HST/ACS im- 
ages of faint galaxies are taken to be perfectly empirical 
representations of reality for the purpose of calibrating 
galaxies dominated by bright HIT regions. 

Furthermore, we have found a good consistency be- 
tween the purely shear-based results (e.g., (-statistic 
based ID reconstruction) and the results based on 
the combined distortion and magnification data (e.g., 
MEM-I- results, B05a results), implying that any shear 
calibration error is not noticeable at the level of our anal- 
ysis, otherwise we would see inconsistency with our mag- 
nification analysis. 

6. SUMMARY AND DISCUSSION 

In this paper, we have derived a projected 2D mass 
map of the well studied lensing cluster A1689 {z = 0.183) 
based on an entropy-regularized maximum-likelihood 
combination of the lens magnification and distortion of 
red background galaxies registered in deep Subaru im- 
ages. The combination of distortion and magnification 
data breaks the mass sheet degeneracy inherent in all 
reconstruction methods based on distortion information 
alone. The method is not restricted to the weak lensing 
regime but applies to the whole area outside the tangen- 
tial critical curve, where non-linearity between the sur- 
face mass-density and the observables extends to a ra- 
dius of a few arcminutes. The strong lensing information 
from ACS observations was also readily incorporated in 
this maximum likelihood approach (§4.4). We also uti- 
lized the distortion measurements to locally downweight 
the intrinsic clustering noise in the magnification mea- 
surements, which otherwise perturbs the depletion signal 
(§4.5). The resulting 2D map showed that the projected 
surface density of A1689 is smoothly varying and sym- 
metric, similar to the distribution of cluster members. 
The 2D mass map is well fitted by the Navarro-Frenk- 
White model, with a continuously steepening profile, but 
the concentration parameter much higher than expected 
for its virial mass (~ 2 x IO^^Mq), according to the clear 
predictions of standard ACDM (Bullock et al. 2001; Neto 
et al. 2007). For consistency we have compared the best- 
fitting NFW parameters obtained from different combi- 
nations of datasets, boundary conditions, and weak lens- 
ing techniques (§5.4.0). We find that all of the recon- 
structions tested here arc consistent with a virial mass 
in the range, Mvir = (1.5 — 2.1) x IO^^Mq, and the com- 
bined ACS and Subaru-2D mass reconstruction yields a 
tight constraint on the concentration parameter, Cyir = 
12.7lJ;9 (^^200 ~ 10), improving upon the statistical accu- 
racy of our earlier ID analysis (B05a). Very good agree- 
ment is found between the present full 2D reconstruction 
and the earlier ID reconstruction (B05a), both of which 
are based on the same distortion and magnification mea- 
surements, supporting the assumption of quasi-circular 
symmetry in the projected mass distribution. We have 
also explored potential sources of systematic error on 
the concentration parameter, suc;h as the uncertainties 
in background redshift distribution, selection criteria for 
the red background sample, and strong lensing model- 
ing. Taking into account all of the systematic errors as 
well as the statistical uncertainty, our constraint on the 
concentration is Cyir = 12.7 ± l(stat.) ± 2. 8 (systematic). 



For clusters well fitted by an NFW profile, the derived 
virial mass of a cluster and the conc;entration parameter 
can be used to find the Einstein radius 9e, through the 
simple relationship 1 = Knfw(^e) (see Appendix B). For 
A1689 with = (2.1 ±0.2) x 10^"^ and c^i, = 12.7^0.9 
(only statistical errors quoted), this yields an Einstein 
radius of 45 ± 6 arcsec at 2^ = 1, or 53 ± 7 arcsec at 
Zs = 3, in very good agreement with the mean estimated 
radius of ~ 50", based on the locations of the multi- 
ple images (Table 2 of B05b). The strong lensing mass 
model of Limousin et al. (2007), based on the multiple 
images identified by B05b, properly reproduces the ob- 
served Einstein radius of 45" at = 1, consistent with 
the strong lensing mass model of B05b (see §5.5.0). In 
contrast, an Einstein radius of 24 ± 11 arcsec at = 1 
is implied by the NFW fit to this cluster by Limousin et 
al. (2007), to independent weak lensing data from CFHT 
(cvir ~ 9.6, Mvir ~ 1.5 X lO^^M©; see Table 5, Figures 15 
and 16). This discrepancy may be attributed to a degree 
of contamination of the sample of galaxies used to define 
the lensed background, which, as pointed out in B05a, 
can drag down the weak lensing signal if accidentally in- 
cluded in the background sample, and preferentially so 
at small radius (see Figure 16) where the ratio of clus- 
ter members compared with background is much higher, 
resulting in a shallow profile and hence a lower con- 
centration fit. 

Taking into account the systematic errors (§5.5), com- 
bined ACS and Subaru constraints allow a shallower (but 
still steeper than theoretically expected) mass profile 
with Cvir = 9 — 10, similar to the values found in Halkola 
et al. (2006) and Limousin et al. (2007). This, however, 
does not simply mean that the discrepancy between dif- 
ferent lensing studies is fully solved: When the NFW 
model is normalized to reproduce the observed Einstein 
radius (45" at Zg = 1), then this concentration woiild 
indicate a large virial mass of Mvir = (3 — 3.3) x W^^AIq 
(see Figures 14 and 15), which however is considerably 
higher than the mass estimates derived from the X-ray 
observations (Mvir ~ IO^'^^Mq in Andersson & Madejski 
2004, Mvir ~ 1.4 X IO^^Mq in Lemzc et al. 2008). The 
discrepancy between the present results and the results 
by Halkola et al. (2006) can be explained by the rela- 
tive weights in the least fitting, assigned differently to 
the ACS- and Subaru-based measurements, indicating 
slight deviation of the observed profile from the NFW 
predictions (see discussion in §5.4.0; also see discussion 
in Medezinski et al. 2007). This is seen at the innermost 
radii r ^ 40kpc//i (Figure 13), where the ACS data indi- 
cate a shallower profile. Consequently, when one prefers 
the innermost region to fit the data, this could lead to 
a lower concentration (cvir = 9 — 10) as favored by the 
shallower slope in the innermost region, and to a lower 
value for the Einstein radius when the data at outer criti- 
cal radiiis are less weighted. Nonetheless, our best-fitting 
NFW model provides a good approximation to our data 
over the radii we have considered, r — [10~^, 2] Mpc/h. 

B05a demonstrate that dilution of the lensing signal 
is certainly the cause of the very low concentration fit 
(cvir ~ 4.5, Table 5) obtained by Bardeau et al. (2005), 
due to the inclusion of relatively blue cluster members 
in the definition of the background sample of the same 
CFHT weak lensing data as Limousin et al. (2007), and 
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for which the equivalent Einstein radius is only ^ 4" (at 
Zs = 1, see Table 5). On the other hand, the dilution 
of the lensing signal caused by cluster members can be 
used to derive the proportion of galaxies statistically be- 
longing to the cluster by comparing the undiluted red 
background distortion signal with the radial distortion 
profile of color-magnitude space occupied by the cluster 
members, but including inevitable background galaxies 
falling in the same space (Medezinski et al. 2007). This 
technique allows the light profile of the cluster to be de- 
termined in a way which is independent of the number 
density fluctuations in the background population, which 
otherwise limit the calculation of the cluster light profiles 
and luminosity functions from counts of member galax- 
ies. The resulting light profile can be compared with 
the mass profile to examine the radial behavior of M/L 
(Medezinski et al. 2007). 

A recent joint X-ray and lensing analysis of A1689 
by Lemze et al. (2008) also produces very similar con- 
centration, Cvir = 12.2^1'^, and virial mass, Mvir ~ 
(1.4 ± 0.4) X IO^^Mq, to that derived here in our analy- 
sis (Table 5). This is derived from a model- independent 
approach to the X-ray emission profile and the projected 
lensing mass profile, assuming hydrostatic equilibrium, 
utilizing the mass profile derived in B05a. Interestingly, 
the observed temperature profile falls short of the pre- 
dicted temperature profile derived from the joint fit and 
this may imply that the gas distribution is clumpy on 
small scales, in the form of a higher density cold gas phase 
(Lemze et al. 2008). Moreover, this anomaly seems to be 
consistent with conclusions regarding gas substructure in 
the recent detailed hydrodynamical simulations of clus- 
ter gas by Kawahara et al. (2007), implying that other 
similar detailed lensing- X-ray studies should also show a 
similar temperature discrepancy. 

Great progress continues to be made in the detailed 
predictions of ACDM, particularly on cluster scales 
where gas cooling is not a worry. Recently the whole Mil- 
lennium survey (Springel et al. 2005) has been converted 
into the observer's frame following the full geodesies 
through the volume to simulate the effect of structure on 
the light received by an observer (Hilbert et al. 2007). 
This work has shown that although in general clusters 
form in overdense regions, the material associated with 
a given cluster in the form of extended groups and fila- 
ments outside the virial radius of the cluster is of rela- 
tively low mass contrast compared to the projected mass 
due to the cluster itself, and therefore lensing based pro- 
jected masses of clusters are not overestimated by more 
that a few percent (Hilbert et al. 2007). 

This simulation has also been used to better define the 
relationship between the concentration parameter and 
the virial mass of halos, over the full range of mass from 
galaxies up to the most massive cluster-sized halos, in the 
context of standard ACDM (Neto et al 2007). A clear 
prediction has emerged that most massive halos gener- 
ated in these simulations have the lower concentration 
(cvir ~ 5), and the cause of this is in part attributed 
to the generally later collapse of the more massive ha- 
los reflecting the lower mean density of the universe. 
For example, at Zvir = 0.183, for the standard choice 
of cosmological parameters the criterion for virialization 
is p(< Tvir) ~ 115/9crit(-Zvir) ~ 277p(2;vir)- One possibil- 



ity to achieve earlier formation of massive clusters is to 
allow deviation from Gaussianity of the primordial den- 
sity fluctuations in the early universe (e.g., Grossi et al. 
2006; Sadeh, Rephaeli, & Silk 2007). 

A degree of triaxiality is inevitable for collisionless 
gravitationally collapsed structures. Discussion of the 
likely effect of triaxiality on the measurements of lens- 
ing properties has been examined analytically (Oguri et 
al. 2005; Sereno 2007; Corless & King 2007) and in nu- 
merical investigations (Jing & Suto 2002; Hennawi et al. 
2007). A bias in favor of prolate structure pointed to the 
observer is unavoidable at some level, as this orientation 
boosts the projected surface mass density and hence the 
lensing signal. This effect has been evaluated in the con- 
text of the CDM model and serves as a guide to the likely 
degree of bias which may affect lensing work. Hennawi et 
al. (2007) conducted a detailed study of the properties of 
lensing cluster population identified in ACDM cosmologi- 
cal A/'-body simulations. The level of bias in terms of the 
concentration parameter derived from 2D lensing mea- 
surements was explicitly estimated and found to amount 
to an ^ 34% increase, which results from a combination 
of two effects, namely the orientation bias (~ 19%) due 
to halo triaxiality and the selection effect towards higher 
3D concentrations (~ 18%). The level of correction for 
the orientation bias is also derived from semi-analytical 
representation of simulated CDM triaxial halos by Oguri 
et al. (2005). The anomalously high concentrations of 
Cvir ^ 13, such as found for A1689, CL0024-M654 (Kneib 
et al. 2003), and MS2137-23 (Gavazzi et al. 2003), ap- 
pear inconsistent with the distribution of concentrations 
found in detailed simulations of Hennawi et al. (2007), 
which predict that only < 2% of lensing clusters should 
have such high concentrations. It is also unlikely that the 
baryonic component in these massive clusters increases 
the concentrations over the ACDM prediction for dark 
matter halos (see discussion in Hennawi et al. 2007, 
B05b, and Broadhurst & Barkanna 2008). 

A chance projection of foreground/background struc- 
ture along the line-of-sight can potentially infiuences pro- 
jected lensing observations, boosting the surface mass 
density locally and hence the concentration. The ACS 
strong lensing analysis of B05b revealed the secondary 
mass clump in the central region of A1689 associated 
with a small clump of galaxies. The existence of this 
subclump has been suggested by earlier observations, 
such as the spectroscopic study of Tcague et al. (1990) 
and Czoske (2004), and the X-ray study of Andersson 
& Madejski (2004). No obvious substructure is visible 
in a large spectroscopic sample of 525 cluster members 
identified in Czoske (2004), in contrast to earlier work 
of Tcague et al. (1990) based on a smaller sample of 
176 identified cluster members. Recently, this secondary 
mass clump has also been directly detected by the weak 
lensing flexion analyses by Leonard et al. (2007) and 
Okura et al. (2008) based on the ACS and Subaru data, 
respectively. However, the detailed ACS strong lensing 
analyses based on multiply imaged background galaxies 
showed that the mass contribution of the secondary mass 
clump is only a small fraction of the main cluster com- 
ponent (Figures 21 and 22 of B05b), implying a lower 
M/L for the subclump than for the main cluster, which is 
tightly constrained by the geometrical positions of sets of 
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multiply lensed images, or the location of critical curves. 

In the near future, the question of the effect of triaxi- 
ality on lensing based cluster mass profiles may be exam- 
ined empirically. For example, the total X-ray luminosity 
or the total lensing based mass of a cluster should not de- 
pend on the orientation with respect to the line of sight, 
whereas the concentration parameter and the Einstein 
radius will be affected and hence expanded lensing stud- 
ies could in principle reveal whether relaxed clusters of 
fixed mass or fixed X-ray luminosity tend to have consis- 
tent lensing based concentrations, or instead a broader 
distribution may be uncovered, with a mean concentra- 
tion smaller than derived for A1689, indicating triaxiality 
produces a significant bias. Current indications based 
on several massive clusters favor the NFW profile, but 
with consistent concentrations (Medezinski et al. 2007, 
in preparation), similar in value to A1689, imderscoring 
the tension between detailed lensing based mass profiles 
and the predictions of standard ACDM. 
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Table 1. Parameters used in the 2D MEM reconstruction 



0.5 96.2 31 X 27 1012 0.700 



''MEM model parameter. 
''Baycsian value of a. 

"^This includes the ACS-constraincd central pixel. 
"^Number of usable measurements. 

''ACS constraint on the central k pixel assuming {D^s/Ds) = 0.693. 



Table 2. 2D MEM reconstructions with different datasets and 

BOUNDARY CONDITIONS 



Method 


Dataset 


Boundary Conditions ^ 




ad 


NDF ° 


2 f 


RMS 8 


S/N h 


2D MEM-S 


shear 


MEM + ACS 


710 


91.2 


515 


448 


0.081 


14.4 


2D MEM 


shear + magbias 


MEM 


1012 


97.3 


785 


1014 


0.076 


18.2 


2D MEM+ 


shear + magbias 


MEM + ACS 


1012 


96.2 


784 


1065 


0.077 


19.4 



Note. — The following three sets of combinations of datasets and boundary conditions are considered: (i) 2D MEM+ method using shear and 
magnification data with the ACS constraint on the central pixel, (ii) 2D MEM method using shear and magnification data without the central ACS 
constraint, and (iii) 2D MEM-S method using shear data with the central ACS constraint. 
''Dataset used for weak lensing analysis. 

''With or without the ACS constraint on the central pixel in the strong lensing region. 
"^Number of usable measurements. 
^Bayesian value of a (m = 0.5). 

"Classical number of degrees of freedom, NDF = N^a-ta, ~ ■'Vgood (see equation [31]). 
'^Minimum functional value of x^- 

S Average rms error of k in the observed region (30' X 24'). 

''Detection significance of the convergence signal defined by equation (38). 



Table 3. Summary of the Methods for Mass Profile Reconstructions 



Method ^ Datasct ^ Boundary Conditions (^mim'^max) ^ ^bin " 



(Jc-statistic 


tangential shear 


Kj, = 4 X 10-3t 


(l',16') 
(l',18') 
(l',18') 
(IMS') 


9 


2D MEM-S 


shear 


MEM 4- ACS 


10 


2D MEM 


shear + magbias 


MEM 


10 


2D MEM-I- 


shear + magbias 


MEM 4- ACS 


10 


Subaru ID 


tangential shear 4 


- magbias ACS 


(l',18') 


10 



''Weak lensing mass reconstruction method. All methods apply to the non-linear but subcritical regime. 
''Datasct used for weak lensing analysis. 

"^With or without the ACS constraint on the central pixel in the strong lensing region. 
^Lowcr and upper radial limits. 

"Number of radial bins in the range of (9min, 9max). 

^This employs an outer boundary condition on the mean convergence Kf, within 16' < < 19'. The mean background level iif, is calculated 
to be Kb = 4 X 10-3 using the ACS4-Subaru-lD best-fit NFW model by B05a. 

ID maximum likelihood analysis by B05a based on the joint measurements of weak lensing distortion and depletion profiles. 
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Table 4. Summary of the best-fitting NFW parameters for Subaru 
weak lensing observations 



Method 


Data 


ACS 


ER 






y2 . /doff 




taiiigeiitiail sneEir 


ID 


9+ 






1 CI +0.27 
^•°^-0.24 




c n /Q /'Q\ 
O.U/ o t^o^ 


Q 

ou.y 


tangential shear 


ID 


9+ 




yes 


1 t;Q+0.24 
1.59_(, 22 


— 2.5 


o 1 /o /'o^ 
y. 1/ y 


44. D 


^c"Statistic 


ID 


K 






1 /I Q 4- n 07 
i.4o It U.Z / 


97 q+2.7 
^'"^-19.3 


o /v /''7^ 
5.2// (0 




Cc-statistlc 


ID 


K 


— 


yes 


1 CI +0.25 


16.5l« 


11.4/8 (8) 


44.8 


Cc-statistic 


ID 


K 


yes 




^ r>i+0.24 
-0.20 


13.7lij 


12.3/19 (19) 


45.3 


2D MEM-S 


2D 


K 








14 1+10.3 


361/834 (406) 


38.2 


2D MEM-S 


2D 


K 


yes 




1 75+0.17 


14.6+1? 


369/846 (418) 


45.1 


2D MEM 


2D 


K 






^•°"-0.17 


14.9lii2^ 


287/834 (478) 


42.7 


2D MEM 


2D 


K 


yes 




1 O1+0.21 
1-^1-0.14 


14.3+1-2 


294/846 (490) 


45.1 


2D MEM-I- 


2D 


K 






-^•'"-0.17 


15.5li:^ 


323/836 (423) 


51.0 


2D MEM-I- 


2D 


K 




yes 


1 93+''-22 

l.»0_Q jg 


14.0l|5 


387/837 (424) 


46.4 


2D MEM+ 


2D 


K 


yes 




2.10 ±0.17 


12 7+1-° 


327/848 (435) 


45.3 


Subaru ID 


ID 


K 








< 30 


5.36/8 (6) 


66.9 


SubaJTu ID* 


ID 


K 


yes 




1.93 ± 0.20 


13.7ll;t 


13.3/20 (18) 


45.4 



Note. — A flat prior of Cvir 5; 30 is adopted in the model fitting. 
'''Mass reconstruction method (see Table 3). 

'^Fitting with or without the ACS-dcrivcd inner mass profile (§5.4.0). 

"^Fitting with or without the strong Icnsing Einstein-radius (ER) constraint, 0^ = 45" at = 1 (§5.4.0). 
"^Virial mass and la error in units of lO^^M©. 
''Concentration parameter, = r^ij-frs, and la error. 

'Values in parentheses refer to effective degrees of freedom excluding upper limit bins with k < 0. 
^Einstein radius in units of arcscc for a fiducial source at = 1, defined as 1 = k{Oe). 
Taken from B05a. A fiducial source redshift of = 1 is assumed in B05a. 



Table 5. Comparison between best-fitting NFW parameters for 
A1689 FROM different observations and methods 



Reference 


Method 


Afvir ^ 


r ■ 


Af200 


£200" 




Remarks 


King, Clowe, & Schneider 2002 
Clowe 2003 
Bardeau et al. 2005 
B05b 


WL 
WL 
WL 
SL 


1.0 
1.3 

1 79+0.8 

3.7 


6.1 
9.9 


0.84 
1.1 

1 41+0.63 

3.2 


4.8 
7.9 

n c + 0.5 

•^•^-o.a 

0-0-1.6 


11 
22 
3.5 
44 


ESO/MPG 
ESO/MPG 
CFH12K 
ACS 


B05a 

Halkola et al. 2006 
Halkola et al. 2006 
Bardeau et al. 2007 


SL-I-WL 
SL 

SL-hWL 
WL 


1.93 ± 0.20 
3.55 ± 0.4 
2.58 ± 0.2 
2.35 ± 0.4 


13.7t\t 
7.6 ±0.6 
9.6 ±0.6 
5.5 ±1.0 


1.72 ± 0.19 
3.05 ± 0.3 
2.25 ± 0.2 
1.97 ±0.3 


6.0 ±0.5 
7.6 ±0.5 
4.28 ± 0.8 


45 
37 
39 
12 


ACS + Subaru {k profile) 
ACS 

ACS + Subaru {g+ profile) 
CFH12K 


Limousin et al. 2007 


WL 


1 -il+o-s 

^•^^-0.2 


9.6 ±2.0 


1.32 ± 0.2 


7.6 ±1.6 


24 


CFH12K 


Lemze et al. 2008* 


SL+WL+X 


2.23 ± 0.6+ 


12.2+?-9 


1.98 


9.7 ±0.8 


45 


ACS + Subaru + Chandra 


This work 


WL 


1.97 ±0.20 


i3.4t-;l 


1.76 ±0.20 


10.7+^-5 


45 


Subaru (k map) 


This work 


SL+WL 


2.10 ±0.17 


12.71-^ 


1.86 ±0.16 


10 1+°'* 


45 


ACS + Subaru {k map) 



Note. — A similar table of best-fitting NFW parameters is found in Comerford & Natarajan 2007 (Table 1) which also include the results for 
other clusters as well as A1689. 
''Analysis method. 

''Virial mass Afvir a-nd Icr error in units of IO^^Mq. 
"^Virial concentration cVir = ^'vir/^'s and la error. 
'^M2oo and la error in units of lO^^M©. 
"Specific concentration C200 = ^"200 /I's and Icr error. 

'^Einstein radius in units of arcsec for a fiducial source at Zs = 1, defined as 1 = k(0e)- 

Based on the Chandra X-ray data and the projected mass profile from the joint ACS and Subaru-ID analysis by B05a. Hydrostatic 
equilibrium assumed. 

^Note that the outermost radius point of Lemze et al. (2008) is at 1.5Mpc/?i, which is smaller than the virial radius of A1689, rvir ~ 
2Mpc//i. As compared to the NFW-based prediction in the table, their model-independent reconstruction of the total mass density p{r) 
yields a virial mass of M^ir = (1-4 ± 0.4) x 10^^ Mq, assuming an extrapolation index of -3. 
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Table 6. Sources of systematic error and their effects on the de- 
termination OF THE halo concentration PARAMETER (Vir- 



Source of error 


Uncertainty range 


Fractional error in Cvir 


Lower color limit for the red sample 


[0.1,0.5] 


4% 


Source redshift 


Zs,D = [0.7, 1] 


10% 


Mask area correction 


^mask = 3 it 1 


7% 


Clustering noise rejection 


I'clust = 4 ± 1 


1.5% 


Slope of unlensed number counts^ 


s = 0.22 ± 0.03 


3.5% 


Inner boundary condition 


K,{Zs = 1) = 0.7810^5 


10% 


Strong lensing modelingtt 


B05b, ER 


10% 


Entropy prior 


m = 0.5 ± 0.4 


10% 



Note. — The stcnialic errors arc presented in fraction of Cvir = 12.7 derived from a joint fit to the ACS-based inner k profile of B05b and the 

Subaru-based k map reconstructed with the MEM+ method. 

Effective source redshift 2s,d equivalent to the mean distance ratio (D^g/Dg) defined by equation (18). 

= dlog N (m) / dm. 

ttEinstein-radius constraint {6e = 45" for Zs = 1) used instead of the ACS inner mass profile of B05b. 
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APPENDIX 

DISCRETIZED ESTIMATOR FOR THE LENSING CONVERGENCE 

In this Appendix, wc aim to derive an expression for the discrete convergence profile using the weak lensing aperture 
densitometry Cc(^) given by equation (46). In the continuous hmit, the averaged convergence R{9) and the convergence 
k{6) are related by 

«W = ^ fdine'e'^Kie'), (ai) 

^(0)-±h^ (A2) 
''^^>-2e^ d\n9 ■ ^^^^ 

For a given set of annular radii 0^ {m = 1,2,..., A^), discretized estimators can be written in the following way; 

2 

K„ = K(e„) = -^^ Aln^,^'/^(^0, (A3) 

™ 1=1 

Ki = K{ei)^alKi+i-a[Ri {l = l,2,...,N-l), (A4) 

where 

-[ = ^('^)\4 = ^('-^]\ (A5) 



2Aln6ii \ eiJ 2Aln( 



with Aln^j = {9i+i — 0i)/9i and 0i being the area-weighted center of the Ith annulus defined by 9i and ^i+i; in the 
continuous limit, we have 

Je, 

_2 9f + ef_^_^ + eiOi+i 

-3 ^^T%^ • ^^^^ 

The technique of the aperture densitometry (Fahlman et al. 1994; Clowe et al. 2000) allows us to measure the 

azimuthally averaged convergence k{9) up to an additive constant k^j corresponding to the mean convergence in the 
outer background annulus with inner and outer radii of 6'inn and 6'out) respectively (§5.3): 

R{9) = Cc{9) + Rb- (A7) 

Substituting equation (A7) into equation (A4) yields the desired expression as 

«(^/) = 4Cc(^i+i) - a[U0i) + (4 - a[)Rb. (A8) 

THE NFW LENS MODEL 
The NFW universal density profile has a two-parameter functional form as 

{r/rs){l + r/rsY 

where Ps is a characteristic inner density, and Vg is a characteristic inner radius. The virial properties are related 
In stead of using and r^, we introduce for an NFW halo the virial mass Mvir and the concentration parameter, 
Cvir = Tvir/rg, defined as the ratio of the virial radius ryir to the scale radius. The virial mass and virial radius are 
related through the following equation: 

47r 

^vir = —p{z^i,)KiKi,, (B2) 

where Avir is the mean overdensity with respect to the mean cosmic density p{z^iv) at the virialization epoch z^i^, 
predicted by the dissipationless spherical tophat collapse model (Peeebles 1980; Eke, Cole, & Prenk 1996). We assume 
the cluster redshift is equal to the cluster virial redshift Zyir. We use the following fitting formula in a flat 3-space 
with cosmological constant (see Oguri, Taruya, & Suto 2001): 

Avir = 187r2(l + 0.4093a;°if 52^, (B3) 

where Wyir = '^/^mizviv) — 1- 
The inner density Ps can be then expressed in terms of other virial properties of the NFW halo: 

Ps = p{ZV.r)^^J— , '''' ^ T. (B4) 

3 ln(l + Cvir) - f;vir/(l + Cvir) 

Hence, for a given cosmological model and a halo virial redshift, we can specify the NFW model with the halo virial 
mass Mvir and the halo concentration parameter Cvir- 



PNFW(»') = \ ^2 (^1) 
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X (pixel) X (pixel) 



Fig. 1. — The quadrupolc PSF anisotropy field as measured from stellar cUipticities before and after the PSF anisotropy correction. The 
left panel shows the raw ellipticity field of stellar objects, and the right panel shows the residual cUipticity field after the PSF anisotropy 
correction. The orientation of the sticks indicates the position angle of the major axis of stellar ellipticity, whereas the length is proportional 
to the modulus of stellar ellipticity. A stick with the length of 5% ellipticity is indicated in the top right of the right panel. 
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Fig. 2. — Stellar ellipticity distributions before and after the PSF anisotropy correction. The left panel shows the raw ellipticity 
components (e*,e2) of stellar objects, and the right panel shows the residual ellipticity components (5e*,5e2) after the PSF anisotropy 
correction. 



For an NFW profile, it is useful to decompose the convergence k{9) and the averaged convergence k{9) as 



H'NFw{x) = -f{x), 



(B5) 
(B6) 



where b = ipsrsfSaitizd, Zs) is the dimensionless scaling convergence, x = 6/{rs/Dd) is the dimensionless angular 
radius, and f{x) and ^(a;) are dimensionless functions. We have analytic expressions for f{x) and ^(a;) as (Bartelmann 
1996): 



fix) 




71=^5 arctanh^ T+f , 



arctan 



rr + l 



(x < 1), 

{x = 1), 
[x > 1), 



g{x)=]ii(^ 



+ < 



7T=f arctanhy' ]^ {x<l), 
1 {x = 1), 

2 0T>^^0-|^ 1 



= arctan (a; > 1). 



The tangential shear 7+,nfw(^) is then evaluated by 

7-|-,NFw(^) = '«NFw(^) — «NFw(^)- 



(B7) 



(B8) 



(B9) 
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Fig. 3. — Averaged shear correction factor, (Pg), as a function of object size, Vg. The horizontal error bar represents the size of the bin 
(Arg = 1 pixel), and the vertical error bar represents the rms scatter in the smoothed scalar correction factor, (Pg). 
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Fig. 4. — Top panel: mean distortion strength averaged over a wide radial range of 1' < 9 < 18', done separately for the blue and red 
galaxy samples. No area weighting is used here to enhance the effect of dilution in the central cluster region. Shown are the measurements 
of the tangential component (3+) with open squares, and those of the 45deg-rotated component (gx)- On the right {red), the square 
symbols show that g+ drops rapidly when the bluer limit of the entire red sample is decreased below a color indicated by the vertical 
dashed line which lies -(-0.22 mag redward of the cluster sequence. This sharp decline marks the point at which the red sample encroaches 
on the E/SO sequence of the cluster. For galaxies with colors bluer than the cluster seqeuence, cluster members are present along with 
background galaxies. Consequently, the mean lensing strengh of the blue sample, as shown on the left {blue), is systematically lower than 
that of the red sample. Bottom panel: the respective numbers of galaxies as a function of color limit, contained in the range 1' < 6* < 18' 
in the red (right) and the blue (left) samples. 



For a given source redshift Zg, the Einstein radius is then readily calculated by 1 = knfw(^); or rnore explicitly, 
using equation (B8) we have 

0| = belg[ej,/es). (bio) 

where Og = rg/Dd — rvir/{cvh-Dd) is the angular size of the NFW scale radius. This equation for 9e can be solved 
numerically, for example, by the Newton-Raphson method. 
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Fig. 6. — Distribution of the lensing magnification bias n/no measured from a red-galaxy sample in the background of A1689, smoothed 
with a Gaussian with FWHM = 2' for visualization purposes. The shaded circle indicates the FWHM of the Gaussian. 
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Fig. 7. — Left panel: The projected mass distribution k of A1689 on a grid of 31 X 27 pixels reconstructed from an entropy-regularized, 
maximum-likelihood combination of Subaru shape distortion and number-count depletion data of red background galaxies. ACS strong 
lensing constraints are used to determine the convergence value at the central pixel that falls in the strong lensing region. The pixel width 
is 1.4 arcmin. The observing region is limited to the central 30' X 24', covering a projected area of 3.9 X 3.1 Mpc^/h'^, while the field size 
for the reconstruction is 43.'4 X 37f8. Overlayed up on the image is the reconstructed spin-2 gravitational shear field. A stick with the 
length of 5% shear is indicated in the top right corner. The north is to the top, and the east is to the left. Right panel: distribution of the 
reconstruction errors for the convergence ft calculated using the maximum entropy method. The uncertainties are given by the square root 
of the diagonal part of the covariance matrix. 



Fig. 8. — Contours {thick) of the dimensionless surface mass density k smoothed with a Gaussian of FWHM = If 4', superposed on the 
V + i' pseudo-color image of A1689. The image size is 30' X 25', covering a projected area of 3.9 X 3.2 Mpc2/fe2 at z = 0.183. The lowest 
contour and the contour interval are 0.05. North is to the top and east is to the left. 
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Fig. 9. — Comparison of the reconstructed lensing fields and the cluster galaxy distributions in A1689. Top: reconstructed convergence 
K {left) and magnification bias fj,^-^''~^ (right) fields. Bottom: observed i'-band surface luminosity (left) and number (right) density 
distributions of (V — i')-selected member galaxies with i' < 23 ABmag. All images are smoothed with a Gaussian of FWHM = 2'. The 
field size is 30' X 24'. North is to the top, east to the left. The contours show the Gaussian smoothed convergence distribution shown in 
the top-left panel. The lowest contour level and the contour interval are 0.05. For each panel the color scale is linear, and ranges from 
0%-100% of the peak value. 
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Fig. 10. — Mass profiles of A1689 obtained by a radial projection of the 2D k map (see Figure 7) reconstructed from an entropy- 
regularized maximum likelihood combination of Subaru distortion and depletion measurements. AH of the profiles are scaled to a fiducial 
source redshift of Zs = 1. The square and triangle symbols represent the results with (2D MEM+) and without (2D MEM) ACS strong 
lensing constrains on the central mass density, respectively. The error bars are highly correlated in the different bins. Without the central 
ACS constraint, central ft-bins at S 4' are slightly underestimated (~ 10%), but the two profiles are overall in good agreement within the 
statistical uncertainties. The solid curve shows the best-fitting NFW profile for the 2D k map (2D MEM-I-) combined with the ACS-derived 
inner mass profile. For comparison an NFW model based on the ACS-|-Subaru ID analysis (1305a) is shown as a dashed curve. 
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Fig. 11. — Comparison of mass profiles from MEM-reconstructed k maps based on different combinations of Subaru datasets and 
boundary conditions. All of the profiles are scaled to a fiducial source redshift of Zj, = 1. The square and triangle symbols represent the 
results from the combined distortion and depletion measurements, with (MEM+) and without (MEM) the ACS constraint on the mean 
surface mass density in the central pixel, respectively. The crosses show the results from the distortion data with the central ACS constraint 
(MEM-S). The model curves are shown for comparison as in Figure 10. The mass profile from the distortion data alone (crosses) shows a 
slight negative dip of k ~ —0.01 at 6' ^ 9 ^ 10 due to spurious boundary effects. 
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Fig. 12. — Shear-based ID mass reconstruction utilizing the (^c-statistic. As an outer boundary condition, the mean background density 
Kf, in the range 16' ^ 6» < 19' is set to K;, = 4 X 10"^ according to the ACS+Subaru-lD best-fit NFW model (B05a). The square symbols 
represent the results with the non-linear corrections. The triangle symbols show the reconstruction with linear approximation. The mass 
profiles are scaled to a fiducial source redshift of Zs = 1- Decorrelated error bars are shown. Downwards-pointing arrows are used where 
the lower error bar drops below zero. Without the non-linear corrections, central bins are underestimated by ~ 15% at maximum. The 
solid curve shows an NFW profile with a high concentration, Cvir = 30, matching well the overall profile obtained with the non-linear 
corrections. For comparison an NFW model based on the combined ACS and Subaru distortion and depletion profiles (B05a) is shown as 
a dashed curve. 
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Fig. 13. — Comparison of model-independent mass profiles of A1689. All of the profiles are scaled to a fiducial source redshift of Zs = 1. 
The filled circles represent the results based on the 2D k map reconstructed from an entropy-regularized maximum-likelihood combination 
of Subaru distortion and depletion data, with the ACS constraint on the mean surface mass density in the central pxiel (MEM-I-). The 
error bars are correlated. The open triangles represent the mass profile from the non-linear fc-statistic measurements based on averaged 
tangential distortion data. Decorrelated error bars are shown. The filled triangles and circles show the results from the ACS strong lensing 
analysis (B05b) and from the Subaru ID weak lensing analysis based on the combined distortion and depletion profiles (B05a). The solid 
curve shows the best-fitting NFW profile for the MEM-reconstructed 2D k map (MEM-I-). The ID- and 2D-based NFW models from the 
respective combined ACS+Subaru data (B05a) are also shown as solid and dashed curves, respectively. For comparison an NFW profile with 
a low concentration, Cvir = 5, normalized to the observed Einstein radius (6e = 45"), is shown as a dotted curve. The low concentration 
model (dotted) predicted for ACDM clearly overestimates the outer profile constrained by the Subaru weak lensing observations. The mass 
profiles are all in remarkable agreement over the full range of radii up to ^ 2Mpc//i. 




Fig. 14. — Joint constraints on the NFW model parameters, (cvir, JV/vir), derived from gravitational lensing observations of A1689. Left 
panels show the 68%, 95%, and 99.7% confidence levels (Ax^ = 2.3, 6.17, and 11.8) in the (c^ir, Mvir)-plane for the 2D ft map reconstructed 
from Subaru weak lensing observations. Right panels show the same confidence levels but for the joint ACS+Subaru-2D NFW fitting, 
incorporating the inner mass profile (10kpc//i ^ r ^ 180kpc/h) constrained by ACS strong-lensing observations by B05b. The cross in each 
panel shows the best-fit set of the NFW model parameters. The top, middle, and bottom panels correspond to the results based on the 
2D MEM-I-, MEM, and MEM-S reconstructions, respectively (see Table 2). The virial mass Mvir is well constrained by the Subaru data 
alone, while the Subaru constraint on the concentration Cvir is rather weak. The complementary ACS observations, when combined with 
the Subaru observations, significantly narrow down the uncertainties on Cvir, placing stringent constraints on the inner mass profile. In 
each panel the observed constraints on the Einstein radius (^e — 45"at Zs = 1) are shown as a dashed curve. 




Fig. 15. — Joint constraints on the NFW model parameters, (cvir, A^vir) obtained from the Subaru tangential shear profile of A1689 
(see Figure 1 of B05a). The cross shows the best-fitting set of the NFW parameters, and the contours show the 68%, 95%, and 99.7% 
confidence levels (Ax^ = 2.3, 6.17, and 11.8) in the (cvir, Mvir)-plane. The observed constraints on the Einstein radius, 6e — 45"at Zs = 1, 
are shown as a dashed curve. The dotted curve shows the Cvir — Afvir relation for 6^ = 24" at Zs = 1. The triangle symbol shows the best-fit 
set of (cvir, A^vir) for the combined ACS and Subaru-2D (MEM+) results. The square and circle show the best-fit sets of (cvir, A^vir) from 
the combined strong and weak lensing analysis of Halkola et al. (2006) and the weak lensing analysis of Limousin et al. (2007), respectively. 
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Fig. 16. — Tangential distortion profile g+{0) {square, upper panel) from the Subaru weak lensing analysis of the red background sample 
(B05a). The solid curve shows the best-fit NFW profile derived from the joint strong and weak lensing analysis of ACS and Subaru 
observations (this work), incorporating full distortion and magnification information. The Einstein radius constraint (triangle) of 0e = 45" 
(zs = 1), determined from multiply lensed images in ACS observations (B05b), is translated to the corresponding depth of the Subaru red 
background sample using the ACS+Subaru-2D NFW model {solid), and added to the distortion profile (g+ = 1), marking the point of 
maximum distortion. The ACS+Subaru-2D NFW model {solid) fits well with the combined ACS and Subaru distortion information over 
the full range of data, r = [80, 2000]kpc//i, but somewhat overpredicts the outer distortion profile at 9 ^ 9' {r Si 1.2Mpc//i). Also shown 
with the dashed curve is the best-fit NFW profile from the CFHT weak lensing analysis of Limousin et al. (2007), which, in contrast, is in 
good agreement with the Subaru outer profile at 8 ^ A' , but underpredicts significantly the inner distortion profile and hence the Einstein 
radius. The dotted curve shows an NFW profile of Halkola et al. (2006) for a simultaneous fit to their ACS inner mass profile and the 
Subaru distortion profile of B05a shown here, but with a different weighting that prefers the inner strong-lensing based profile where the 
data imply a shallower slope (see Figure 13). The lower panel shows the radial profile of the 45° rotated component gx {9) for the same 
Subaru red background sample (B05a). The X -component of the red galaxy sample is consistent with a null signal at all radii, indicating 
the reliability the Subaru distortion analysis. 



